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ABSTRACT 

We present new Hubble Space Telescope Cosmic Origins Spectrograph far-ultraviolet (far-UV) spec¬ 
troscopy and Keck Echellete optical spectroscopy of 11 ultraluminous infrared galaxies (ULIRGs), a 
rare population of local galaxies experiencing massive gas inflows, extreme starbursts, and promi¬ 
nent outflows. We detect Lya emission from 8 ULIRGs and the companion to IRAS09583+4714. 

In contrast to the P Cygni profiles often seen in galaxy spectra, the Lya profiles exhibit prominent, 
blueshifted emission out to Doppler shifts exceeding -1000 km s _1 in three HH-dominated and two 
AGN-dominated ULIRGs. To better understand the role of resonance scattering in shaping the Lya 
line profiles, we directly compare them to non-resonant emission lines in optical spectra. We find that 
the line wings are already present in the intrinsic nebular spectra, and scattering merely enhances the 
wings relative to the line core. The Lya attenuation (as measured in the COS aperture) ranges from 
that of the far-UV continuum to over 100 times more. A simple radiative transfer model suggests the 
Lya photons escape through cavities which have low column densities of neutral hydrogen and become 
optically thin to the Lyman continuum in the most advanced mergers. We show that the properties 
of the highly blueshifted line wings on the Lya and optical emission-line profiles are consistent with 
emission from clumps of gas condensing out of a fast, hot wind. The luminosity of the Lya emission 
increases non-linearly with the ULIRG bolometric luminosity and represents about 0.1 to 1% of the 
radiative cooling from the hot winds in the HH-dominated ULIRGs. 


1. INTRODUCTION 


Hydrogen recombination in star-forming galaxies pro¬ 
duces prodigious amounts of line radiation including 8 
Lya photons for every Ha photon. The emergent Lya 
and Ha line profiles typically look very different be¬ 
cause galaxies are nearly always optically thick to the 
Lya photons. Even when the neutral gas column be¬ 
comes optically thin to Lyman continuum radiation, i.e., 
N(HI) < 1 /cro ~ 1.59 x 10 1 ' cm -2 where <To represents 
the hydrogen photoionization cross section at the Lyman 
edge, the Lya optical depth at line center remains very 
large, 


r„(Lya)=8.02xlO»(iApU) 

/ N(HI) \ 
\1.59 x 10 17 cm -2 / 


( 1 ) 


using the atomic data in IMortonl (1200311 and a Doppler 
parameter b = \j2kT/ ijih- The Lya photons random 
walk out of the galaxy and undergo diffusion in frequency 
space with each scattering. The increased path length 
makes the line photons more susceptible than the under¬ 
lying continuum to absorption by dust. The luminosity 
and frequency of the emergent Lya photons are there¬ 
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fore sensitive to the dust content, phase str ucture, and 
kinematics of interstellar gas (|Diikstrall2014f) . 

We can better understand these properties of the in¬ 
terstellar medium by studying Lya emission from nearby 
galaxies for which many physical properties have been 
measured, and, as a result, also learn how to best de¬ 
rive physical properties of higlr-redshift galaxies from 
Lya properties when few other diagnostics are possible. 
Given the large number of fac tors influencing the radia - 
tive transfer of Lya photons (iGronke fe Diikstral I2014T ) . 
empirical studies are needed for galaxies sampling the 
full gamut of dust properties and gas kinematics. 

Prior to the Galaxy Evolution Explorer (GALEX) 
spectroscopic survey, limited access to the space ultravi¬ 
olet restricted observations of Lya 1 215.6 7 A to less than 
two d ozen ga laxies at low redshift (IGiavalisco et alJll996t 
iKunth et al.ll!998h . A systematic search of GALEX data 
produced 66 Lya- emit ting galaxies at 0.2 < z < 0.35 
(jDeharveng et all [20 08?! : and a study of nine deep fields 
extended the reds hift range of GALEX Lya emitters 
(ICowie et al.ll2010t ). Multiwavelength follow-up demon¬ 
strated that these low-redshift, Ly q-em i tters, h ave un- 
dergone very recent star formation (ICowie et al.ll201 Th . 

These studies have introduced the concept of a Lya 
escape fraction, f esc (Lya), which compares the flux de¬ 
tected in the Lya line to that predicted from the extinc¬ 
tion corrected Ha flux und er Case B reco mbination con¬ 
ditions at low density. lAtek et ahl (12009H demonstrated 
the systematic decline in / esc (Lya) with increasing dust 
extinction. The attenuation of the resonance line in 
many galaxies exceeds the attenuation of the underly¬ 
ing far-UV continuum. This result suggests resonance 
scattering increases the pathlength of the Lya photons 
and therefore the probability of absorption by dust grains 
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(j Atek et al.1 I2014D . Resonance scattering also explains 
the larger angular size of the gal axies in Ly a relat ive to 
the population of massive stars (lHaves et al.ll2013l) . 

The recent ad ditio n of t he Co smic Origins Spectro¬ 
graph (COS, IGreen et al.l (120121 ) to the Hubble Space 
Telescope (HST) substantially improved ultraviolet sen¬ 
sitivity. This new opportunity has allowed the ex¬ 
ploration of the Lya properties for a wider range of 
galaxy types dHeckman et al.l 1201 ll iWofford et all 120131 : 

iLeitherer et al.l 1201. 'll IJaskot & Oevl 120141 )' Particu¬ 

larly relevant COS observations for understanding Lya 
emission from high-redshift galaxies include the pio¬ 
neering studies_of_Jow-redshift Lyman Break Analogs 
(LBAs. [Heckman et al.1 (120111) 1 and of lower mass galax¬ 
ies wit h extreme specific star formation rates (Green 
Peas. ETaskot fo Oev (120141) b New COS observations of 
the most luminous, dusty starbursts (introduced and dis¬ 
cussed here) and a forthcoming analysis of a larger sam¬ 
ple of Green Pea galaxies (Henry et al., in prep.) build 
on this work. 

The properties of some LBAs and Green Pea galax- 
ies make them candida tes for Lyman continuum leakage. 

I Heckman et al.l (2001) showed that sightlines optically 
thin to the Lyman continuum would not generate com¬ 
pletely black troughs in low-ionization state absorption 
lines. The discovery of significant residual intensity in 
the C II A1334.5 feature of both t he L BAs with the mos t 
dominant central objects (DCOs. lHeckmai]_et_alj j2011JV) 
and a few Green Pea galaxies (iJaskot fc Oev 12014?) 
theref ore generated much interest, and IBorthakur et al.l 
(120141) subsequently measured an escape fraction of a few 
percent from the direct detection of the Lyman contin¬ 
uum from the DCO J0921+4509. The observations of 
both populations also suggested the Lya emission prop¬ 
erties were modified in s ystems wit h esc a ping Lyman 
continuum. In particular, iHeckman et alJ ( 20111 ) noted 
a correlation between the presence of blueshifted Lya 
absorption and a high (inferred) escape for the Lyman 
continuum. 

None of these studies have observed the most lumi¬ 
nous sites of star formation in the local universe, de¬ 
terred perhaps by their enormous infrared excess. Dusty 
starbursts, however, comprise an increasing fraction of 
the c osmic star formation rate (SFR) with lookback 
time dLe Floc’h et al.ll2005l) ; and, somewhat surprisingly, 
sub-mm selected galaxies at high redshift actually show 
Lya in emission more frequently than d o Lyman-break 
selected galaxies (Chapman et al.1 f2005Tl . Observations 
of the Lya emission from local Ultraluminous Infrared 
Galaxies (ULIRGs) are largely missing from the lit¬ 
erature yet would provides important insight into the 
physics of Lya radiative transfer and feedback in intense 
star-forming regions at a much earlier time. 

In this paper, we present new COS observations of 11 
galaxies selected on the basis of their high far-infrared 
luminosity, Lib > 10 12 L 0 , but representing a range of 
stages in the transformation of gas-ric h sp iral gala xies t o 
field ellipticals via major mergers (iGenzel et al.l [2001 1. 
As a class, these objects emit a tiny fraction of their 
luminosity in the far-UV, but their enormous luminos¬ 
ity makes them detectable with COS. We expect sub¬ 
stantial variations among their absorption and emission 
properties in Lya due to both the wi de variation in the 
UV/FIR flux ratio among ULIRGs (Murray et al.ll2007h 


and the evolution of the interstellar medium (ISM) dur¬ 
ing the merger. Because of their warm dust temperatures 
and high luminosities during the latest stages of merg¬ 
ers, a subset of ULIRGs have been recognized as quasars 
partly hidden by a l arge central concentration of dust 
(Sanders et al.1 [1*98 8 alibi ). 

In Section [21 we describe the selection of the ULIRGs, 
the new HST observations, and the supporting optical 
spectroscopy. We characterize the Lya radiation from 
ULIRGs in Section [3] emphasizing the unusual shape of 
the line profile. Our interpretation of the Lya profile is 
enhanced by a quantitative comparison to recombination 
and forbidden lines in optical spectra in Section [I] We 
discuss the physical origin of the broad emission in star- 
burst galaxies in Section [5] and present simple radiative 
transfer modeling of the Lya line profiles in Section [6] 
The results point towards the highly blueshifted emission 
coming from a physically distinct component of the out¬ 
flow related to the hot wind fluid, and we summarize the 
evidence for this results and its implications in Section [7] 

We adopt a cosmology with H 0 = 70 km s -1 Mpc -1 , 
U 0 = 0.3, and Uy = 0.7 throughout. At the median 
galaxy redshift of z = 0.11, one second of arc subtends 
2.0 kpc. 


2. OBSERVATIONS 

We have previ ously used the Eche llete Spectrograph 
and Imager fESI lSheinis et alJ (120021 )1 on the Keck II 
telescope to obtain high- resolution optical spectra of 
many of the 2 Jy U LIRGs (Straus s et al.l ll992l) imaged in 
K and R bands bv lMurphv et al.l ( 19961) . We found out¬ 
flows of low-ionization-state gas via interstellar Na I and 
K I absorption lines that are blueshifted several hu ndred 
km s' 1 with respect to the systemic redshift (Martini 
I2005L [20061) . We also resolved broad wings on many op¬ 
tical emission lines and demonstrated that shocks trav¬ 
eling at several hundred km s -1 excite this componen t 
of the emission dSoto et al.1 l2012t iSoto fc Martinll2012l) . 
To gain further insight into these gas flows, we have ob¬ 
tained new G130M spectra of ULIRGs from the above 
sample with the Cosmic Origins Spectrograph (COS) on 
the Hubble Space Telescope (HST, GO programs 12533 
and 13407). 


2.1. The COS ULIRGs Sample 
2.1.1. UV Selection 

We selected 11 of the UV-brightest ULIRGs at red- 
shifts from 0.084 to 0.151 for the COS observations. Like 
all ULIRGs, our targets have extremely high ratios of far- 
infrared to far-UV flux. Most of the radiation from their 
stars and AGN is absorbed by dust, and the energy re¬ 
radiated as a thermal spectrum with maximum intensity 
between 60 //in and 100 //to. Their bolometric luminosi¬ 
ties are just 0.30 dex brighter than their far-infrared lu¬ 
minosities which we list in TableQ]for reference. Very lit¬ 
tle of their bolometric luminosity comes out in the ultra¬ 
violet. By choosing the UV-brightest ULIRGs, however, 
our sample is biased towards systems with lower than 
average (for ULIRGs) color excess, and this selection ef¬ 
fect does we aken the average strength of int erstellar Na I 
absorption (Martinll2005HChen et al.ll2010D . 

The bolometric luminosities of the COS ULIRGs span 
a relatively narrow range from log(L/jj/ L 0 ) = 12.07 to 
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12.55 b ut are corr e lated with the merger stage. Follow¬ 
ing [viilleu^etjil] (|2002D . we estimate the merger state 
from the separation of the nuclei and prominence of the 
tidal tails as described in the notes to Column 8 of 
Table [l] The most luminous object, IRAS11598-0112, 
has progressed to stage IVb of the merger sequence. 
The lowest luminosity targets (IRAS09583+4714 and 
IRAS08030+5243) are wide binaries still in the pre¬ 
merger stage (Ilia) with two distinct galaxies still recog¬ 
nizable. The far-infrared luminosity of IRAS09583+4714 
is dominated by one of the two galaxies, which we call 
IRAS09583+4714; but it is the companion 14.5 kpc 
away which emits most of the far-UV flux. We pointed 
COS at this companion which we will refer to as 
IRAS09583+4714-B. 

The properties of the COS ULIRG sample are also con¬ 
sistent with the emergence of an AON at late times. 
Most of the COS ULIRGs have the cool far-infrared 
color typical of starburst dominated ULIRGs. How¬ 
ever, the relative flux densities at 25 /xm and 60 fi m, 
which are listed in column 4 of Table [U identify two ob¬ 
jects - IRAS01003-2238 and IRAS 12071-0444 - with the 
higher dust temperatures indicative of a more compact 
source. The upper limit on the color for IRAS11598- 
0112, F v (25fj,m)/F l/ (60/xm) < 0.22, also allows a clas¬ 
sification as a warm ULIRG. These warm ULIRGs are 
less common than cool ULIRGs and mark the tran¬ 
sition from cool ULI RGs to optically selected quasars 
(jSanders et al.l[l988bfl . The three warm ULIRGs in our 
sample show compact morphologies indicative of a late 
stage merger. For example, IRAS01003-2238, is the 
most compact merger in the sample (stage V). Both 
IRAS12071-0444 and IRAS11598-0112 are classified mor¬ 
phologically as (stage IVb) compact mergers. 

2.1.2. The AGN Contribution to LboI 

In recent years, much progress has been made towards 
understanding the AGN contribution to the bolometric 
luminosity of ULI RGs. In particular, the ULIRG evolu¬ 
tion study of lVeilleux et all ( 2009 ), used six independent 
methods to estimate the AGN fraction and found consis¬ 
tent luminosity fractions to within 10-15%. We can use 
these results to estimate the AGN contribution to the 
bolometric luminosity of each COS ULIRG. 

For example, Figure |T| shows the optical spectrum 
of IRAS11598-0112. The narrow forbidden lines and 
the broad Balmer e missi on lines ide ntify t his UL IRG as 
a Seyfert 1 galaxy (iOsterbrock & FerlancU l2Ci5t^ ). The 
presence of strong Fe II emission in the optical spec¬ 
trum de fines this object as a Narrow Line Seyfert 1 
(NLS1) (iLeiglrlv fe Mooreil2004lh The ULIRG evolution 
study showed that AGN provide roughly 75% of Ljn in 
Seyfert 1 ULIRGS, so the luminosity of IRAS11598-0112 
is likely dominated by the AGN. 

The emission-line ratios measured from the optical 
spectra of IRAS01003-2238 and IRAS12071-0444 place 
them in the Seyfert 2 region of the excitation diagrams 
(ISoto V Martini f2(Ti~2lh The Seyfert 2, IRAS01003-2238, 
is the most advanced merger (stage V) in the sample; 
and, since this is consistent with the scenario where the 
AGN emerges towards the end of the merger, we consider 
this object AGN dominated. The emission-line nebula 
in IRAS 12071-0444 is larger than t he typical narrow -line 
region of an AGN, however, and (|Soto et al.l I2012H pre- 





A (A) 


Fig. 1. — Optical spectra of the Seyfert 1 ULIRG IRAS 11598- 
0112 zoomed in on features that identify this object as a narrow¬ 
line Seyfert 1 galaxy. In addition to broad Balmer lines and more 
narrow forbidden lines, the optical spectrum shows strong Fe II 
emission complexes. 

viously attributed the excitation to shocks related to a 
galactic wind rather than an AGN. Still, he high far- 
infrared luminosity of IRAS 12071-0444, its warm color, 
and late (stage IV) merger class ar e all consistent with a 
significant AGN contribution. The IVeilleux et al.1 (|2009f) 
results find that AGN power as much as 50% of the to¬ 
tal luminosity from Seyfert 2 ULIRGs, we infer that the 
AGN luminosity could be as large as the starburst lu¬ 
minosity in both IRAS01003-2238 and IRAS12071-0444. 
For most of the COS ULIRGs, their HH-dominated op¬ 
tical spectra and cool far-infrared color strongly suggest 
the starburt domina tes the bolometric lu minosity. For 
this class of ULIRG, IVeilleux et al.1 (|2009fl estimate that 
AGN contribute no more than 15% to 35% of the bolo¬ 
metric luminosity. 

2.2. HST Imaging & COS Target Acquisition 

For ULIRGs with single-nuclei, we centered the bright¬ 
est UV knot in the 2'/5 diameter primary science aper¬ 
ture (PSA) using the dispersed-light acquisition proce¬ 
dure recommended by the COS Instrument Handbook 
for Cycle 19. In those ULIRGs that have double nuclei, 
however, the UV emission may be dominated by star 
formation in either galaxy or even compact knots in the 
tidal tails (|Surace fe Sanderil2000H . For all but the com¬ 
pact ULIRGs associated with the later merger stages, 
we obtained short exposures with HST through far-UV, 
near-UV, and optical filters and found the coordinates 
for the UV-brightest knots prior to acquiring each target 
with COS. Table [2] lists the parameters of these imaging 
observations. 

The standard pipeline processing for each imaging con¬ 
figuration proved sufficient for target acquisition pur¬ 
poses. To construct images suitable for photometry, how¬ 
ever, we applied a zero-level bias correction to the raw 
images and ran LACosmic to remove cosmic rays before 
combining frames with AstroDrizzle. After these steps, 
significant trails can still be seen on the detector from 
the incomplete transfer of charge toward the amplifier. 
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Figure [2] shows these optical, near-UV, and far-UV im¬ 
ages. The ACS SBC/F125LP and WFC3 UVIS/F225W 
images generally detect the same knots, but many of 
the morphological features present in the optical F625W 
images remain completely undetected in the ultraviolet. 
Even the optical images can be blind to the most ac¬ 
tive regions which are v eiled by thick layers of dust, e.g., 
Ay ~ 50 dGenzel et al.lll998h . 

The far-UV images illustrate the distribution of light in 
the COS Primary Science Aperture (PSA) for our spec¬ 
troscopic observation. The sizes of the ULIRGs range 
from 0"17 to 1//07 FWHM, smaller than the full diam¬ 
eter of the PSA but similar to the unvignetted portion. 
For example, the COS PSA subtends 5.03 kpc at the 
median target redshift; but the diameter of the unvi¬ 
gnetted region is just 1.61 kpc. It is unknown whether 
a low surface-brightness halo of scattered Lya emission 
extends beyond the PSA for some targets. 

Our strategy for computing blind offsets used all 3 im¬ 
ages. The F125LP bandpass provides the best match 
to our spectroscopic bandpass, so we used it to iden¬ 
tify the brightest knot; but the field of view is typi¬ 
cally too small to include a bright point-source. The 
WFC3/UVIS F225W images provide a larger field-of- 
view which typically includes a star or quasar bright 
enough for a near-UV imaging acquisition; and the blind 
offsets were computed from these images. The short 
WFC3/F625W exposures provide easy registration with 
optical images and were used to verify absolute pointing 
and estimate the red leaks in the UV filters. In one case, 
IRAS09583+4714, the stars in the F225W image were 
not bright enough for a blind offset. Following the same 
procedure adopted for ULIRGs without pre-imaging, we 
centered the PSA on a bright knot significantly offset 
from the position of the maximum far-infrared intensity, 
hereafter IRAS09583+4714-B. 

2 . 3 . COS Far-UV Spectroscopy 

We obtained COS spectra in TIME-TAG mode 
through the PSA using the medium resolution far-UV 
grating, G130M. We observed at four focal plane offset 
positions in order to reduce the impact of fixed-pattern 
noise associated with the micro-channel plate. The spec¬ 
tra covered much of the bandpass from 1137 to 1468A. 
The Lya line fell on Detector B in the observation of the 
lowest redshift target, IRAS23365+3604, but was red- 
shifted onto Detector A for all other observations. In 
each case, the redshift of the galaxy cleanly separated 
the Lya profile from geocoronal emission. Table [3] lists 
the central wavelengths (CWLS) of each pointing and 
additional details about the observations. 

Individual frames were processed with version 2.13.6 of 
the CalCOS pipeline. Caicos corrects the data for instru¬ 
mental effects, assigns a vacuum wavelength scale, and 
extracts flux-calibrated spectra. It applies a heliocentric 
correction to the final xld files for each exposure. We 
combined the xld files usi ng soft ware developed by the 
COS GTO team (|Danforth et al.lfeOlOf) . A simple mean 
pixel combination gave better results than the default 
variance weighting for these low S/N data. 

2 . 3 . 1 . Wavelength Calibration 

Emission from geocoronal airglow fills the COS aper¬ 
ture, so the observed wavelengths of these lines are inde¬ 


pendent of target position in the PSA. We measured the 
wavelengths of Lya 1215.67, 01 1302.1685, 01 1304.8576, 
OI 1306.0286, OI] 1355.5977, and OI] 1358.5123 but did 
include NI 1199.967 line because it is a blend of N I 
1199.55, 1200.22, 1200.71. The velocity offsets among 
these lines were always < 40 km s -1 , which we adopt as 
the maximum error in the pipeline dispersion solution. 

Any error in centering the galaxy in the COS aperture 
will generate systematic errors in the wavelength calibra¬ 
tion. To verify the centering of the target in the PSA, we 
identified Milky Way absorption lines in the spectra with 
continuum detections. An offset of the target from the 
center of the aperture will systematically shift the wave¬ 
lengths of these lines. In the spectra IRAS01003-2238, 
IRAS11598-0112, IRAS16487+5447, IRAS00262+4251, 
and IRAS09583+4714, we measured five lines - Si II 
1190.42 Si II 1193.29 Si II 1260.42 Si III 1206.50 C 
II 1334.53; and we measured a subset of these in the 
IRAS23365+3604 and IRAS12071-0444 spectra. The 
continuum S/N ratio is not adequate to measure any of 
these lines in the IRAS10378+1109, IRAS17028+5817, 
and IRAS08030+5243 spectra. Blends with other tran¬ 
sitions are the main source of error in the velocity off- 
sets@ After rejection of the obvious blends, the average 
velocity offset of the Milky Way absorption was within 
the uncertainties of the wavelength calibration for most 
targets. 

We identified three spectra, however, with Galactic 
absorption lines that were Doppler shifted by amounts 
greater than the uncertainty in the wavelength calibra¬ 
tion. We compared these Doppler shifts to the IWakkerl 
(l2004fl map of Galactic high velocity clouds and made 
the following associations. Since the only significant ve¬ 
locity offsets we found have the same sign as the Doppler 
shifts of high-velocity clouds in their direction, we con¬ 
cluded that the UV knots were well centered in the COS 
PSA and did not apply any corrections to the pipeline 
wavelength scale. Hence when comparing the veloci¬ 
ties of UV and optical lines, we consider features within 
±40 km s" 1 to be at the same velocity. 

For completeness, we also note that our sightline 
towards IRAS00262±4251 passes 3.2° northwest of 
M 31, a projected distance of just 44 kpc. In the 
IRAS00262±4251 spectrum, we detect Si III, Si II, and 
S II at approximately -30 0 km s -1 , an offse t consistent 
CGM (see, for example. iLehner et all (120141 1). 

2 . 3 . 2 . Spectral Resolution 

The G130M resolution is excellent when observing 
point sources, R ~ 18,000; but degrades to R ~ 1450 
(about 207 km s -1 ) when the angular size of the target 
exceeds the diameter of the COS PSA. The effective res¬ 
olution of each COS ULIRG spectrum depends on the 
angular size of the UV knots but must lie between these 
limits. 

6 The Si II 1260.42 can blend with S II 1259.52. In the 
IRAS00262+4251 spectrum, for example, the other two lines in 
the S II 1250.58, 1253.81, 1259.52 triplet are detected and un¬ 
blended making the S II 1259 blend with Si II 1260.42 very obvi¬ 
ous. The blend gives the Si II 1260 absorption trough an anoma¬ 
lous blueshift. The CII 1334.53 trough is blended with the weaker 
CII* 1335.71 line in some spectra, and this blend systematically 
increases (redshifts) the velocity offset. This CII blend lies on top 
of the IRAS00262+4251 Lya line in that spectrum. 
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Fig. 2.— HST images of the ULIRGs in the optical (WFC3/F625W), near-UV (WFC3/F225W), and far-UV (ACS/SBC/F125LP) bands 
plus an archival image of IR.AS01003-2238 in WFPC2/F439W. Targets are ordered from the highest to lowest Lyo luminosity. The images 
are registered with axes labeled in arcseconds relative to the coordinates of the COS pointing given in Table [3] north is up and east is the 
the left. The diameter of the COS aperture is 2'/5, but the flux is heavily weighted by the central 0”8 diameter region due to vignetting. 
The galaxies are compact in the ultraviolet bands relativ e to their optica l imag es. The UV knots are typically offset by a few tenths of an 
arcsecond from the peak optical light consistent with the lOoldader et ab (2002 i study. 
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We estimated the actual resolution of our spectra 
two ways. We first measured the widths of the un¬ 
blended Milky Way absorption lines in the 5 spec¬ 
tra with good continuum S/N ratio. Then we esti¬ 
mated the resolution from the size of each source in 
the near-UV imageQ These measurements had four ob¬ 
jects in common, and we found the Galactic absorp¬ 
tion lines always indicated lower resolution by a fac¬ 
tor of two to three relative to the images. We at¬ 
tribute this to the extended nature of the UV emis¬ 
sion which is not captured by the FWHM. Using the 
Galactic lines, the COS G130M resolution for the most 
compact ULIRG, IRAS01003-2238, is 69 km s -1 while 
the stage Ilia merger IRAS 16487+5447 fills the COS 
aperture and degrades the resolution to approximately 
200 km s" 1 . For IRAS16474+3430, IRAS00262+4251, 
and IRAS11598-0112, we measure a spectral resolution of 
150 km s , 154 km s -1 , and 159 km s _1 , respectively. 
Based on the size of IRAS09583+4714-B in the near-UV 
image, we would estimate the resolution R < 1830 or, 
equivalently, FWHM > 164 km s -1 . 

2.3.3. The Ultraviolet Continuum 

When the ultraviolet continuum was detected, we typ¬ 
ically found spectral signatures of young stars. Broad, 
shallow O VI and N V absorption troughs accompanied 
by redshifted emission, for example, indicated the pres¬ 
ence of stellar winds from very massive stars. The ul¬ 
traviolet continuum is solidly detected in six of the COS 
spectra, and only the spectrum of IRAS11598-0112 fails 
to directly indicate a population of ionizing stars. 

Since the UV spectral slope of a young stellar popula¬ 
tion is more sensitive to reddening than to age, we could 
estimate the steller reddening Eb-v,i from the observed 
shape of the UV continuum. We first c orrected each spec- 
trum for Galactic reddening using the lFitznatrickl (|1999T ) 
extinction curve. Next, we assumed an intrinsic UV spec¬ 
tral slope appropriate to a y oung stellar populatio n, i.e., 
/3 ~ —2.3 where F\ oc A' 3 ([Leitlierer et al 1 IT9991) . and 
reddened t his template using a Calzetti extinction law 
(ICalzetti et al.ll2000f) . Using custom software which al¬ 
lowed spectral binning and interactive adjustments to the 
fitting window, we fit the steller reddening Eb~v ,* and 
the UV continuum simultaneously. 

The uncertainty in our estimate of the stellar reddening 
was taken to be the larger of (1) the range resulting from 
placement of the continuum or (2) the range obtained by 
repeating the exercise with /3 = —2.1 and /3 = —2.5. The 
resulting values of the stellar reddening in Table 0] range 
from 0.02 to 0.32. These values are not large compared 
to normal galaxies and demonstrate that the UV light 
emerges through holes in the overall dust distribution. 

The global continuum fits generally seemed reason¬ 
able around individual spectral lines. Where the pseudo¬ 
continuum formed by the blending of lines in the under¬ 
lying galaxy/AGN spectrum was significant, however, we 
fit a linear continuum to bandpasses on either side of the 
spectral line and normalized the spectrum by the fitted 

7 We have measured the profile of the near-UV surface bri ght- 
ness and then used the scaling given bv IHeckman et ahl 1)20111 1 for 
COS observations of the DCO’s, i.e., FWHM of O'.'l to 0'.'5 yields 
a resolution R = 13, 000 to 3000 with the medium resolution grat¬ 
ings. 


continuum. We illustrate the allowed range for the con¬ 
tinuum placement near Lya in Figure [3] We used these 
limits to estimate the uncertainties in our measurements 
of the Lya emission. 

2.4. Echellete Optical Spectroscopy 

We previously obtained optical spectra with ESI in 
echellete configuration on_the Keck II telescope as 
described previously dMartinl l2005f) . The resolution 
through the 1"0 by 20'.'0 ESI slit, R ~ 5000 or 60 km s -1 
FWHM, is comparable to that obtained for the most 
compact ULIRGs with COS G130M. These slits inter- 
se ct the locations of th e COS pointings and are shown 
in ISoto fe Martini (I2012D . We extracted an ESI spectrum 
along a 2'.'5 length of each slit centered at the position 
of the COS aperture. The rows of the spectrum were 
weighted by their distance from the center of this aper¬ 
ture using the vignetting function for the G130M grat- 
ing@ 

2.4.1. Emission-Line Diagnostics 

All the ULIRGs show strong Balmer and forbidden 
line emission in these ESI spectra. We fit these emis¬ 
sion li nes simultaneously using the procedure described 
in ISoto et al.1 (|2012D which includes a stellar Balmer ab¬ 
sorption component. Table [5] gives the reddening esti¬ 
mated from the Balmer decrement. This nebular color 
excess is similar to normal galaxies because it does not 
probe the dust in the regions completely opaque to op¬ 
tical light. The reddening towards these HII regions is 
generally, but not exclusively, larger than our rough esti¬ 
mates of the stellar reddening from the UV spectral slope 
in Section 12.3.31 as would be expected considering that 
the younger ages of the stars emitting the Lyman con¬ 
tinuum limit the distances the stars have traveled from 
t heir birth clo uds. _ 

ISoto fe Martini (120121 1 previoulsy measured well- 
known, emission-line diagnostics from the ESI spec¬ 
tra and placed pairs of emission-line ratios in dia¬ 
grams that distinguish excitation by photoioni zatio n 
from m assive stars from radiative shocks (All en et al.1 
l2008fl . ISoto et, all (|2012f ) found the broad component 
of the line profile is excited exclusively by shocks mov¬ 
ing at speeds of 200-400 km s -1 . In Table 0 we list the 
shock velocity estimated from the line ratios [OIII]/H/3 
and [OI]/Ha at the position of the COS aperture. 

2.4.2. Doppler-shifted Resonance Absorption 

Measurements of the Na I resonance absorption for 
IRAS00262+4251, IRAS08030+5243, IR AS16487+5447 , 
and IRAS23365+3604 were collected from lMartinl (|2005T) . 
Table [5] shows these equivalent widths and new mea¬ 
surements made in a similar manner for the other 
COS ULIRGs and IRAS09583+4714-B. The ESI spec¬ 
tra of IRAS 10378+1109 and IRAS08030+5243 show the 
strongest Na I absorption at 6.7 and 7.4 A, respec¬ 
tively. The Na I doublet in IRAS09583+4714-B and 

8 COS ISR-2010-09 provides throughput measurements in only 
the XD (cross dispersion) direction for G130M. The Instrument 
Science Report compares the vignetting function for G160M in 
the XD and spectral directions and shows that they are nearly 
identical; hence, we have tacitly assumed the G130M vignetting is 
also quite similar the XD and spectral directions. 
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IRAS17028+5817 is not strong enough to definitively 
claim the detection of a non-stellar component. The 
interstellar Na I resonance absorption, when detected 
however, is blueshifted roughly 200 km s -1 at maximum 
depth, and the blue wing becomes undetectable around 
-700 km s^ 1 as it approaches the continuum level and 
blends with a He I emission line. An exception is 
IRAS08030+5243, the__only examp le of redshifted Na I 
absorption shown in iMartinl (1200511 . 

2.4.3. Aperture Corrections for Ha Luminosity 

Most applications of the optical spectral measurements 
rely on line ratios rather than absolute fluxes, so the ge¬ 
ometrical difference between the circular COS aperture 
(2"5 diameter) and the rectangular ESI aperture (1"0 by 
2"5) is not a concern, especially since the COS through¬ 
put falls rapidly at radii beyond 0"4 from the center. 

The geometrical distinction between the rectangular 
ESI slit and circular COS PSA needs to be taken into 
account, however, when we compare the s trengths of Ha 
and Lya. We used the R band images dMurphv et al.1 
1199611 to model the Ha aperture correction because the 
scale lengths of ga lactic disk s in the R band and in Ha 
are very similar dKennicutti 11989 ). Specifically, we vi¬ 
gnetted the R image using the COS function described 
above, and then measured the R band flux in the 2"5 
diameter COS aperture relative to that in a 1"0 by 2"5 
rectangular ESI aperture. To estimate the Ha flux that 
would be measured through a 2"5 diameter aperture, we 
multipled the exctinction-corrected Ha fluxes by these 
aperture corrections which we also list in column 4 of 
Table [5] 

The spatial extent of Lya halos can also introduce sys¬ 
tematic errors in esti mates o f f e ^(Lya). For example, 
as demonstrated in lHaves et all (12013), resonance scat¬ 
tering enlarges galaxies in Lya by a factor of 1 to 3.6 
over their size in Ha; hence, our aperture may include a 
larger fraction of the total Ha emission than of the to¬ 
tal Lya emission. This type of bias is unlikely for the 
COS ULIRG sample, however, because the LARS data 
empirically show that a low dust content is a necessary 
requirement for Lya halos__to exte nd well beyond the re¬ 
gion of Ha emission (lHaves et al.ll2013l ). 

3. COS SPECTROSCOPY OF Lya 

Figure [3] shows the COS G130M spectra around Lya. 
We detect emission from eight of the ten ULIRGs. 
Considering the large optical depth of the Lya tran¬ 
sition, see Eqn. [I] and extreme concentration of dust 
in ULIRGs, it is remarkable that any Lya photons es¬ 
cape from ULIRGs. It is perhaps easier to understand 
the absence of Lya emission from IRAS0803+5234 and 
IRAS17028+5817; the upper limits on their Lya fluxes 
are far below those detected. 

To better appreciate the surprising shape of the 
ULIRG Lya profiles, we show the spectrum of the non- 
ULIRG, IRAS09583+4714-B, in Figure® The Lya emis¬ 
sion is redshifted roughly 250 km s -1 relative to the op¬ 
tical emission lines, and the resonance absorption in sev¬ 
eral Si II transitions and Si III is blueshifted 100 km s -1 
or more. Blueshifted Lya absorption is visible after 
smoothing the spectrum. P Cygni line profiles with red¬ 
shifted emission and blueshifted absorption as seen in 
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Fig. 3.— Rest-frame Lya emission in COS G130M spectra of 
ULIRGs. Airglow lines are denoted by ©, and the rest wavelength 
of Lyais indicated by the thick tick marks at the top and bottom. 
Spectra are ordered from the highest to lowest line luminosity (de¬ 
noted as a logarithm in ergs s~ 1 ). The Doppler shift of the blue 
wing, Vb 20 is given in km s —1 . The fitted continuum and its un¬ 
certainty range are shown by the horizontal solid and dashed lines, 
respectively; note the strong pseudo-continuum from metal lines in 
the spectrum of the mini-BAL IRAS11598-0112. With the excep¬ 
tion of IRAS23365+3604, the blue line wing reaches higher velocity 
as the ULIRG luminosity rises. 


this spectrum arise from r esonance scattering in a n ex¬ 
panding, neutral medium (IVerhamme et all 1200 8). Six 
of the 13 Ly a emitters from the Hct-selected galaxy sam¬ 
ple of IWofford et al.1 (120131) and all four of the infrared 
luminous galaxies (10 11 < Lir/ L© < 10 11 ' 7 ) observed 



























































bv lLeitherer et all (120131) show these P Cygni Lya pro¬ 
files; and they have been modeled with outflowing shells 
of neutral gas expanding at velocities < 200 km s -1 . 

The dominance of blueshifted emission among the 
ULIRG Lya profiles indicates different physical condi¬ 
tions than those which produce the P Cygni profiles. The 
profiles are unlike the local Ha-selected galaxies without 
P Cygni profiles; they display Lya absorption troughs 
with damping wings fit by H I column d ensities in the 
range of 1.0 x 10 20 to 2.3 x 10 21 cm -2 (|Wofford et al.1 
12013 ). Even among UV-selected galaxies a t z = 2 — 3 
where 30% show multiple emission peaks, iKulas et all 
(120121 ) find the redshifted peaks are normally stronger 
than the blueshifted emission component. Only the line 
profiles of IRAS16474+3430 and IRAS16487+5447 show 
most of the net emission redward of the systemic velocity, 
and only IRAS16474+3430 may (depending on the con¬ 
tinuum placement) have a weak, blueshifted absorption 
component. 

We have ordered the ULIRG spectra by decreasing 
Lya luminosity in Figure[3] The most luminous ULIRG, 
IRAS11598-0112, is a quasar on account of its dominant 
AGN and bolonretric luminosity; blending of its broad 
Lya profile with emission lines from N V, Si II, and C III 
produces line wings much broader than on any other COS 
ULIRG. If a G160M spectrum of this object is obtained 
in the future, the C IV emission could be used as a tem¬ 
plate pro file to debl end this complex as shown, for ex¬ 
ample, in lLeighlv et all ( 2007 ). For now, the line flux we 
measure from integrating this broad emission complex 
provides an upper limit on the Lya luminosity due to 
the contamination from other lines. 

Blueshifted emission wings are prominent on the 
Lya profiles of IRAS01003-2238, IRAS 12071-0444, 
IRAS10378+1109, and IRAS23365+3604. These wings 
extend roughly -1000 km s _1 from the systemic veloc¬ 
ity. The process that generates these line wings does 
not appear to be uniquely associated with an active nu¬ 
cleus. Although we argued in Section 12.1.11 that the 
starburst and AGN make roughly equal contributions 
to the bolometric luminosities of IRAS01003-2238 and 
IRAS12071-0444, starbursts power IRAS 10378+1109 
and IRAS23365+3604. Radiative transfer modeling has 
previously linked blueshifted Ly a emission to either ra¬ 
dial infal l dDiikstra et~al . 2006a, b j) or low columns of neu- 
tral gas (jVerhamme et al.ll2014l iBehrens et al.ll2014D . 

3.1. Escape Fraction of Lya Photons 

To gain insight into the role of radiative transfer effects 
in shaping the Lya emission, we compare the strengths 
of the Lya and Ha emission. Scattering off neutral hy¬ 
drogen atoms increases the pathlength of Lya photons 
through the ISM. The Lya photons will be more atten- 
tuated than the Ha radiation due to a combination of 
the longer distance traveled due to scattering and the 
frequency dependence of interstellar reddening. We will 
distinguish between these processes using the attenua¬ 
tion of the far-UV continuum. 

The intrinsic ratio F(Lya )/ F{ Ha) depends only on the 
electron density and temperature for Case B conditions, 
i.e., the ionized region is optically thick to photons just 
above the hydrogen edge {IP = 13.60 eV). In the low 
density limit, we have n e « n e , C rit = 1.55 x 10 4 cm" 3 
and collisions can be neglected. The intrinsic Lya flux is 
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Fig. 4.— Continuum normalized spectra of the galaxy 
IRAS09583+4714-B. The COS G130M spectrum of the 
IRAS09583+4714-B galaxy show a P Cygni Lya profile consistent 
with an outflowing wind or shell. The UV absorption lines from 
low-ionization-state gas are blueshifted relative to the systemic 
velocity set by the optical emission lines. The deepest part of 
the absorption trough is shifted roughly 100 km s —1 whereas the 
blueshifted wing extends to no more than 300 to 400 km s —1 . 
Other lines detected, but not shown, include Si II 1260.42 which 
is blended with S II 1259.52, S II 1253.81, Ar I 1066.66, N II 
1083.99, and multiple Fe II lines. The optical spectrum shows 
no absorption in the Na I doublet indicating sodium is mostly 
ionized. 


then 8.1 times the intrinsic Ha flux (lDraineH2011 ). 

Galaxies typically have Lya to Ha flux ratios much less 
than the Case B ratio. It is therefore common to describe 
the suppression of Lya photons relative to the intrinsic 
Lya flux indicated by the Ha emission. We define the 
escape fraction of Lya photons as 


fesc (Lya) = 


F(Lya) 


8.1 x F(Ha) c 


( 2 ) 
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We calculate the intrinsic Ha flux from the observed 
Ha flux and the nebul ar reddening u sing the relation 
appropriate for a lCardelli et al.l (119891) reddening curve, 
namely 


F cor (Ha) = F obs ( Ha) x 10 1012£b - v 


The color excess for the ionized gas is 


Eb-v 


1 


0.9692 


log 


F(Ha)/F(H£) 

2^86 


(3) 

(4) 


and was computed from the same reddening curve. 
We t a ke the intrinsic Balrner decrement from 
lOsterbrock fc Ferlandl (120061 ) for Case B conditions, an 
electron density of n e = 100 cm -3 , and a temperature 
of T e = 10 4 K. 

The Lya escape fraction computed this way ranges 
from 1-3% among the starburst-dominated spectra with 
Lya emission. For our two non-detections, we find 
/ esc (Lya) < 0.5%. Figure [5] shows that these escape 
fractions decline with increasing color excess. Among 
the nine Lya detections, the Spearman rank order corre¬ 
lation coefiicent is rg = 0.42, and the correlation deviates 
from the null hypothesis of no correlation by just 1.2cr. 
These escape f raction s fall within the range measured for 
local galaxies (lAtek et al.ll201 4) as indicated in Figure 0 

To assess the importance of resonance scattering, we 
compare the attenuation of Lya to that of the contin¬ 
uum emission near Lya. The continuum attenuation is 
shown as a function of the color excess by the dashed line 
in Figure 0 From this comparison we see that the Lya 
attenuation in three of the ULIRGs - IRAS1207-0444, 
IRAS10378+1109, and IRAS23365+3604 - is similar to 
the continuum attenuation. In spectra of the other non- 
AGN ULIRGs and IRAS09583+4715-B, the Lya escape 
fractions is greatly suppressed relative to the continuum 
due to scattering. The upper limits on the Lya flux for 
the two undetected ULIRGs imply the resonance line is 
attenuated over a factor of 100 more than the far-UV con¬ 
tinuum. We conclude that the Lya optical depth varies 
widely among the COS ULIRGs, scattering has the least 
affect on the Lya emission emerging from IRAS1207- 
0444, IRAS23365+3604, and IRAS10378+1109. 


3.1.1. High f esc (Lya) from a Mini-BAL Quasar 

Our measurements allow unusually high values of 
f esc (Lya) for IRAS11598-0112, but the accuracy of this 
estimate merits additional discussion because a quasar 
powers this object. The uniqueness of this object’s 
far-UV spectrum among the COS ULIRGs might have 
been anticipated from the st rong optical Fe II emission 
lines s hown in Figure ITl which iLioari et all (I1993D : ILioaril 
(1994) have previously associated with the emergence of 
nascent broad absorption line quasars in ULIRGs. Our 
detection of broad, high-ionization-state lines in Figure[6] 
solidly identifies this object as a mini-BAL (broad ab¬ 
sorption line) quasar. 

The assumptions underlying the intrinsic strength of 
Lya in Equation [2] may not hold for IRAS11598-0112. 
For example, the intrinsic Case B ratio of F(Lya)/F(Ha) 
will be larger than 8.1 as n e approaches the critical den¬ 
sity; collisions promote electrons from the 2s state to 
2 p where they can decay by emission of a Lya photon. 
In some quasar spectra, the intrinsic ratio is driven in 
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Fig. 5.— Escape fraction of Lya photons vs.(nebular) color ex¬ 
cess. The inverted triangles are upper limits. The dashed line 
illustrates the suppression of th e continuum near Lya from dust as 
given by a lCardelli et all 11989 ) reddening curve and the color ex¬ 
cess (measured from the Balrner decrement). When the Lya escape 
fraction of a galaxy lies far below this curve, we attribute the ex¬ 
tra Lya attenuation to the longer path length caused by resonance 
scattering. The ULIRGs do not stand out from other samp les of 
Lya emitting galaxies shown here for comparison. See Section [3. 1.11 
for further discussion of the Seyfert I galaxy, IRAS11598-0112. 

the oppposite direction (i.e., towards enhanced Balrner 
emission re lative to L ya) by the rmalization of the en¬ 
ergy levels (jKwan fc Kroliklll979D . Thermalization only 
becomes important at very high densities however. 

The far-UV spectrum illustrated in Figure 0 helps 
us place constraints on the electron density. The 
IRAS11598-0112 outflow is only detected in lines from 
high-ionization-state gas. We detect blueshifted absorp¬ 
tion in the O VI 1032, 1038, S IVA1062, and N V 1239, 
1243 transitions; and the widths of these lines place 
IRAS11598-0112 in the mini-BAL (broad absorption 
line) quasar category. To directly measure n e , we need to 
measure the abso rption from a collisionally excited leve l 
of one of the ions (iBorguet et al.ll20l3l : lDunn et al.ll2012T) . 
Our spectrum covers S IV* A1073 at comparable S/N ra¬ 
tio to the resonance transition at A1062 A. No absorption 
is detected from the excited line, so the electron density 
lies well below the critical density of 3.16 x 10 4 cm -3 . To 
summarize, the density is too low for thermalization to 
determine the level populations of hydrogen but is likely 
high enough for the intrinsic Case B ratio to be signifi¬ 
cantly larger than 8.1. 

By assuming the low density limit, we have underes¬ 
timated the value of the denominator in Eqn. [2] As ex¬ 
plained above, however, we have overestimated the Lya 
flux in the spectrum (i.e., the numerator in Eqn. 0) by 
including blended metal lines in our Lya flux. It fol¬ 
lows that we obtain no quantitative limit on / esc (Lya) 
for IRAS11598-0112. Our point is that the fiducial value 
stands out compared to the other ULIRGs, and we at¬ 
tribute this to the dominant AGN in IRAS11598-0112. 
We speculate that this object is viewed through an ion¬ 
ized cone (with neutral walls) pointed in our direction. 
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Fig. 6. — COS G130M and optical M/1 spectra of the narrow-line, 
Seyfert 1 ULIRG IRAS11598-0112. (a) The Lya and Ly/3 profiles 
show strong absorption near the systemic velocity. The hydro¬ 
gen Balmer lines have broad wings reaching nearly 2000 km s~ 1 . 
While we expect such wings on the Lya profile, the even broader 
wings visible in the top panel are believed to be heavily blended 
with broad emission from metal lines including Si II. See also the 
[O III] profile in the Appendix, (b) The COS G130M spectrum 
shows strong absorption from high-ionization gas. No absorption 
is detected from low-ionization-state gas or Si III. The absorption 
from high-ionization state gas is blueshifted up to 500 km s —1 . The 
equivalent widths of O VI A1032 and N V A1239 are 1.43 A and 
1.29 A, respectively. 


In radiative transfer calculations, the escape fraction of 
Lyct photons from a compact source in an anisotropic gas 
distributions can be boosted by factors of 2 — 4 relative 
to the isotrop i c case when viewed through low-N(HI) 
holes (jBehrens et al.ll2014fl . 


3.1.2. Correlation of f esc ( Lya) with the Width of the Lya 

Profiles 

To facilitate direct comparison of the linewidth with 
f esc ( Lya), we need a quantitative characterization of 
broad wings and/or line asymmetries. Describing the 
broad wings and asymmetries of the Lya profiles ro¬ 
bustly, however, provides a challenge. The full width of 
the profile at half maximum intensity, for example, is not 
uniquely defined for the double peaked profiles, nor does 
it flag line profiles with broad wings or significant asym¬ 
metry. In Appendix A, we define velocity markers Vb 
and Vr as a continuous function of the fractional profile 
area. For any percentile of the profile area, the full width 
of the interpercentile velocity (IPV) is well defined. The 
widths at the 0.20, 0.50, and 0.80 percentiles are marked 
on each Lya profile in Appendix Figures a-i to illustrate 
this concept. 

The escape fraction increases as the width of the Lya 
emission line increases. This correlation is both more sig¬ 
nificant and stronger than that for f esc { Lya) and color 
excess. The correlation with line width reflects an un¬ 
derlying correlation of / esc (Lya) with the Doppler shift 
of the blue wing. This correlation is stronger than the 
one with total linewidth, and / esc (Lya) does not show a 
correlation with the velocity of the red wing. Figure [7] 
compares / esc (Lya) to the Doppler shift of the blue wing. 

This measured correlation is in the opposite direction 
to that expected for a static slab. For a static slab, the 
width of the resonance line would increase with the num¬ 
ber of scatterings. Since the longer pathlengths would 
decrease / esc (Lya), we would expect the Lya escape frac¬ 
tion to decline as the line width increased. 

Perhaps the simplest intrepretation of this correlation 
is that / esc (Lya) depends strongly on the gas kinemat¬ 
ics, not just the column density of gas and dust. How¬ 
ever, previous empirical evidence for this effect ha s been 
limited to studies of a few galaxies dKunth et all 119981 
lAtek et all 120081 : iWofford et al.l 120131 ) . Among the HII- 
dominated ULIRGs, we find / esc (Lya) is correlated more 
strongly with emission-line width than with color excess. 

3.1.3. The Lya Equivalent Width 

Often the Ha spectrum is not accessible for Lya emit¬ 
ters. In this situation, the Lya equivalent width may be 
a useful indicator of the attenuation of Lya photons. The 
equivalent width measures the strength of the emission 
line relative to the continuum^ The ULIRGs as a class 
show moderate equivalent widths similar to the values 
measured previously for the local IUE-selected galaxies 
plotted Figured] The median IF (Lya) is 11 A among 
the COS ULIRGs, and the strongest equivalent width 
ULIRGs include the AGN. 

For continuous star formation, the UV continuum (like 
the Ha emission) scales linearly with the SFR. In con¬ 
trast to the F(Lya)/F(Ha) ratio, however, the strength 
of Lya relative to the UV continuum is quite sensitive to 
the age of the stellar population because the UV contin¬ 
uum can be produced by stars of insufficient mass to gen¬ 
erate a substantial ionizing continuum. Considering the 

9 We define the Lya equivalent width by an integral over the 
line such that W(Lya) = f(F\ c ( A) — F\(\))/F\ c {X)d\, where 
F\ c is the fitted continuum. It represents the continuum bandpass 
(measured in A) with the same flux as the emission line. 
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Fig. 7. — The Lya escape fraction increases with the blueshift of 
the Lya emission wing, defined here as the velocity beyond which 
10% of the Lya emission line flux is measured. The Spearman rank 
order correlation coefficient is r s = 0.83, and the null hypothesis 
is rejected at 2.4<r. The correlation strength is driven by the large 
fesc (Lya) values inferred for the ULIRGs with the highest AGN 

fraction. _ 

sensitivity of W( Lya) to stellar age dCowie et al.11201 1.) 
then, the expected strength of the correlation between 
these Lya measurements is not clear a priori. 

Figure [5] illustrates how |lU(Lya)| grows with 
f esc (Lya). Even among the starburst-dominated 
ULIRGs, there is a nearly one-to-one correspondence, 
and we conclude that / esc (Lya) and fU(Lya) indicate 
a consistent picture of moderate Lya escape fractions 
from starburst dominated ULIRGs. In the opposite pic¬ 
ture where stellar age dictates W{ Lya), the ULIRGs fur¬ 
thest along the merger sequence would in fact be ex¬ 
pected to have the oldest stellar populations and there¬ 
fore the weakest IU(Lya). In contrast, we find the op¬ 
posite. The most advanced mergers tend to be AGN 
dominated ULIRGs, and these objects show the highest 
Lya equivalent widths. 

3.2. Caveats 

In this section, we have shown that Lya escape frac¬ 
tions from ULIRGs are comparable to those in more typ¬ 
ical star-forming galaxies. We explain this surprising re¬ 
sult by the fact that the COS spectroscopy is blind to 
the most active and heavily obscured regions in ULIRGs. 
Unlike normal star clusters and starburst galaxies which 
show a strong correlation between the far-in f rare d ex¬ 
cess and the UV spectral slope dMeurer et al.ll 1995 1. the 
ULIRGs have a very blue UV spectral slope for their 
infrared excess. No UV radiation would be detected 
from ULIRGs if the dust distribution could be approx¬ 
imated by the foreground screen model which works so 
well for local galaxies. To illustrate this point, we used 
the GALEX magnitudes and far-infrared luminosities 
from Table HI and es timated the infrared excess following 
iMeurer et all (119951) . We find the typical infrared excess 
value in the COS ULIRGs is IRX = F fir /F 220 ~ 10. 
We have only measured the escape fraction relative to 


Fig. 8. — Correlation of f esc (Lya) with equivalent width of Ly« 
emission. The ULIRGs with the highest escape fractions also have 
large equivalent widths. The correlation coeffient is r s = —0.82 
for the nine detections, and the null hypothesis is rejected at the 
2.3 ct level. The lower limits on |VF(Lya)| arise from the absence of 
continuum detections. The downward pointing triangles illustrate 
the upper limits on the two ULIRGs with no Lya detection; these 
galaxies were arbitrarily placed at zero equivalent width since the 
continuum is not detected either. 

the regions transparent to optical light. 

To estimate the escape fraction is in an absolute sense, 
we compare the SFR from unobscured regions, as mea¬ 
sured from the Ha luminosity (not corrected for internal 
extinction), to the SFR indicated by the far-infrared lu¬ 
minosity. We estimate that unobscured regions account 
for only 0.3 to 2 percent of the total star formation rate 
among the COS ULIRGs, and the measured escape frac¬ 
tions refer only to these unobscured regions. The abso¬ 
lute escape fractions are far smaller than those estimated 
from UV and optical lines. 

The COS PSA subtends 5.0 kpc at the median redshift 
of the sample, and measurements of Lya flux from the 
COS spectra will underestimate the total Lya luminosity. 
A more specific concern relevant for the / esc (Lya) results 
is that resonance scattering can produce large Lya ha¬ 
los around galaxies relative to the a ngular extent of the 
Ha emission. Extrapolation of the lHaves et all (|2013t ) 
results for the LARS galaxies (which found the largest 
differences in spatial extent among galaxies with very low 
dust abundance) would suggest little spatial bias for the 
COS ULIRGs sample. Considering the unusual condi¬ 
tions in the ISM of ULIRGs, it would be quite interesting 
to directly image the spatial extent of their Lya emission 
and test whether they follow the trends for halo emission 
defined by more normal galaxies. Since the far-UV im¬ 
ages have peak sensitivity near the wavelength of the Lya 
emission, this may be possible using using sophisticated 
spectral fitting. techniques as outlined, for example, by 
lHaves et all (j2009P . 

3.3. Implications 

Our discovery of blue wings on the Lya line profiles 
raises the question of whether ULIRGs might have gas 
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flows with bulk velocities of roughly ~ 1000 km s -1 . The 
large Doppler shift of the Lya photons does not necessar¬ 
ily require fast outflows. The frequency diffusion which 
accompanies repeated scatterings can produce a large 
Doppler boost. For example, a resonance photon scat¬ 
tered back and forth repeatedly between opposite sides 
of an expanding shell could get a Doppler boost that 
greatly exceeds the shell velocity. Distinguishing bulk 
motion from frequency diffusion is the first critical step 
towards understanding the conditions that produce the 
line wings. 

From our analysis of f esc (Lya) in Section 13.11 we 
concluded that the Lya photons from IRAS12071-0444, 
IRAS23365+3604, and IRAS10378+1109 do not undergo 
many scatterings because their attenuation is similar to 
that of the UV continuum; hence, they probably do not 
receive their large Doppler boost from frequency diffu¬ 
sion. We also find that / esc (Lya) is most strongly cor¬ 
related with the Doppler shift of the blue wing on the 
Lya profile, the Lya line width, and the Lya equivalent 
width. 

We explore these concepts further in Section 0] We 
use the optical line profiles of the ULIRGs as a template 
for the intrinsic profile of the emitted Lya photons, i.e., 
prior to their passage through the galaxy, and discuss the 
implications for the Lya radiative transfer. 

4. COMPARISON OF H LYMAN AND OPTICAL LINE 
PROFILES 

The optical emission lines in the Keck II/ESI spectra 
of the ULIRGs provide a view limited to the optically 
transparent holes. This perspective is, however, closely 
related to the regions visible in the ultraviolet spectra; 
many common features can be seen in the images pre¬ 
sented in Figure [5] for example. In Figure [H]a,b, we show 
the line profiles of the H Balmer lines and [O III] A5007 
forbidden line on the same velocity scale as the COS 
G130 Lya emission. We begin this section by comparing 
their line cores and then demonstrate the relationship be¬ 
tween the highly blueshifted gas seen in the Keck II/ESI 
and COS G130M spectra. 

4.1. Profile Core 

In Figure [HI the velocity of the Lya absorption 
typically coincides with the velocity of the maxi¬ 
mum optical emission which is near the systemic ve¬ 
locity. This absorption in the profile core is rela¬ 
tively narrow and most prominent in the HH-dominated 
ULIRGs (IRAS0026+4251, IRAS10378+1109, and 

IRAS16474+3430) as well as the three systems with 
prominent AGN (IRAS11598-0112, IRAS01003-2238, 
and IRAS 12071-04444). The Lya absorption follows our 
expectation that neutral hydrogen in these galaxies will 
scatter Lya photons away from the systemic velocity. 

One unexpected feature of these line cores is the 
residual intensity at the systemic redshift. In the 
IRAS16487+5447 spectrum, the central absorption 
troughs lie to either side of the peak Ha intensity. In 
IRAS11598-0112, we see emission in Lya in the broad 
absorption trough near the velocity of the peak optical 
emission. We also note that many of the spectra are 
not completely black in the absorption troughs. In a 
few cases, this residual intensity near line center may re¬ 
sult from the line-spread function smearing out narrower 


lines which are black at line center. The troughs that are 
completely black (IRAS11598-0112, IRAS09583+4714, 
IRAS16474+3430, and IRAS16487+5447) appear signif¬ 
icantly broader however, and our highest resolution spec¬ 
trum, IRAS01003-2238, shows a narrow trough which is 
not black. The residual intensity requires some sightlines 
with very low column densities of neutral hydrogen, a re¬ 
sult we will refine quantitatively in Section [6] 

Scattering of resonance photons out of the line core 
will broaden the Lya line profile. Following absorption 
near the resonance, the Lya photons random walk in 
frequency space until they are either absorbed by a dust 
grain or are Doppler shifted to a frequency where they 
can escape. We therefore expect the Lya lines to be 
somewhat broader than the optical lines. 

To investigate the relationship between the observed 
and intrinsic Lya profiles, we use a non-resonant opti¬ 
cal emission-line to approximate the intrinsic Lya profile 
shape. We have normalized the Lya and [O III] line 
profiles by their respective areas and overplotted them 
in Appendix A. To make a quantitative comparison, we 
use the technique based on profile area introduced in Sec- 
tion !3. 1.21 and further described in Appendix [XI We mark 
the interpercentile velocity at p = 0.80 and take this as a 
measurement of the width of the core of the line profile. 
In the plot of IPV as a function of profile area, compar¬ 
ison of the Lya and [O III] line widths at p = 0.80 illus¬ 
trates the systematically broader width of Lya. Further 
inspection of thsese plots shows that IPV(Lya) exceeds 
IPV([0 III]) at all values of the fractional profile area p. 
It follows that the Lya profiles are systematically broader 
than the optical line profiles regardless of exactly where 
the width is measured. The only ULIRG with an [O III] 
profile broader than Lya is IRAS16474+3430 which has 
been broadened by orbital motion; the 2D version of the 
ESI spectrum shows a second nucleus offset l'/84 from 
the center of the COS pointing. We therefore conclude 
that diffusion of resonance photons in frequency space 
typically explains the broadening of the cores of the Lya 
profiles relative to the non-resonant lines. 

4.2. Line Wings 

In the optical lines shown in Figure [9] the hydro¬ 
gen Balmer lines typically show a broad, blueshifted 
line wing. The median blueshift of the fitted wings is 
250 km s -1 dSoto et al.ll2012H . but this emission compo¬ 
nent often extends to 1000 km s -1 due to its large veloc¬ 
ity dispersion. This high velocity gas emits strongly in 
optical forbidden lines of [O III], [S II], [N II], and [O I]. 
We focus the discussion on the comparison between the 
[O III] 5007 line and Lya because the former lies in a 
spectral region free of blends with other transitions. In 
contrast, if we use Ha or H/3, then the deblending of 
[N II] 6548, 84 or, respectively, stellar H/3 absorption in¬ 
troduces additional noise into the analysis. 

The only measurable consequence of this choice is that 
our template for the optical emission line profile has 
slightly stronger wings relative to the line core than it 
would had we used a Balmer line. The shape of these 
optical profiles differs slightly due to a change in exci¬ 
tation mechanism. Line ratio diagrams show that the 
broad component is not excited by the ionizing contin¬ 
uum from massive stars. Since its spatial extent is also 
typically larger than the broad line region of an AGN, 
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ijSoto et al.l I2012T) have previously fitted the diagnostic 
line ratios measured for the broad emission component 
with models for radiative shocks. 

Figure HU] compares the extent of the blue wings of the 
[O III] and Lya emission. For purposes of illustration, 
we chose the interpercentile velocity markers at p = 20% 
as a compromise between being away from the line core 
yet not too sensitive to continuum placement. The ve¬ 
locity of the blue Lya wing is correlated with that of 
the [OIII] wing. We also see that the measurements of 
F/320 (Lya) and Vb2o([OIII]) scatter around the locus of 
equality. Even though the cores of the line profiles are 
relatively broader in Lya, the blue wings are not any 
more extended than the [O III] line wings at the same 
fractional profile area. This result is noteworthy because 
the line profiles of the non-resonant line are not affected 
by scattering. The good agreement suggests the intrinsic 
profile of the Lya emission contains blue line wings. 

The Doppler shifts of the blue wings therefore arise at 
least in part from bulk motion along the line-of-sight and 
cannot be attributed solely to frequency upscattering in a 
static medium. Since the same flows are detected in the 
optical emission lines, the Doppler shift of the optical 
lines provides as good an indicator of / esc (Lya) as did 
Lya in Figure [7] Figure ITT! shows this correlation of the 
Lya escape fraction with the velocity of the blue wing 
on the [O III] profile. 

Inspection of the area normalized profiles in the Ap¬ 
pendix enables a closer comparison of the blue wings in 
Lya and [O III]. The more similar the line profiles, the 
less radiative transfer shapes the Lya profile. We draw 
particular attention to subfigures d (IRAS10378+1109), 
f (IRAS12071-0444), and i (IRAS23365+3604). Between 
velocities of roughly -200 km s _1 and -1000 km s -1 , the 
shapes of their Lya profiles are nearly identical to that 
of the optical line profile. We therefore conclude that 
radiative transfer effects play little or no role in shap¬ 
ing the blue wings in these ULIRGs. This result fur¬ 
ther strengthens our conclusion that the bulk velocity of 
the Lya emitting gas must reach at least 1000 km s _1 in 
these objects. 

In stark contrast to the blue wings, the velocities of the 
red wings on the Lya profiles are not correlated with the 
optical line profiles. Inspection of the area-normalized 
profiles in Appendix Figures a-i shows that the red wings 
on the Lya line profiles extend to larger Doppler shifts 
than do the [O III] profiles. The greater extent of the red 
wing in the resonant line makes the linewidths systemat¬ 
ically larger for Lya. The red wings on Lya also weaken 
the correlation between the Lya and [O III] linewidths FI 
Due to the presence of the red wings, the Lya profiles are 
actually less symmetric than the optical line profiles by 
some measures. 

For example, in the Appendix, we define a line asym¬ 
metry measurement across the line profile (i.e., computed 
at each interpercentile velocity). The asymmetry com¬ 
pares the relative strength of the emission redward and 
blueward of the median velocity and is positive for a net 
blueshift. The panels for IRAS12071-0444 illustrate this 
concept for a Lya profile with two maxima. Near the 

10 The linewdiths are not quite as strongly correlated as are the 
velocities of the blue wings. At IPV 20 for example, we measure 
rs = 0.85 with the null hypothesis rejected at the 2Acr level. 


core (p> 0.4 in this spectrum), the red velocity marker is 
further from the median velocity than is the blue veloc¬ 
ity marker at the same IPV; however, in the line wings 
( p< 0.3) the blueshifted profile has the larger Doppler 
shift. Since the optical lines do not have the redshifted 
components in their line profiles, their net asymmetry 
parameter is positive and biased towards the blue. The 
redshifted component of the Lya emission is typically not 
present in the optical line profiles and therefore provides 
another profile metric that successful models for the Lya 
emission must explain. 

Indeed, the main conclusion from our comparison to 
non-resonant line profiles is that resonance scattering 
does not create the high velocity wings. The absence 
of red wings on the optical line profiles suggested this 
emission arises from an outflow near the starburst re¬ 
gion. In this scenario, we see emission from blueshifted 
gas on the near side, but dust in the ULIR G b locks the 
redshifted emission from the far side dSoto et al.l I2012T ). 
Measurements of line profiles at wavelengths that prop¬ 
agate through the dust generally support this interpre¬ 
tation. The CO 1-0 line profile of IRAS23365+3604, for 
example, exhi bits br oad (albeit very weak) symmetric 
line wings (jCicone et all 12014 ). It therefore seems likely 
that the blue wings observed on the Lya line profiles 
are directly emitted by outflowing gas, and dust between 
us and the far side of this outflow attenuates the red 
wings on the intrinsic Lya profile. The broad Lya wings 
must be emitted by gas moving at velocities of 500 to 
1000 km s -1 . 

4.3. Comparison to Absorption-Line Kinematics 

The fast winds indicated by the line wings present 
a paradox because their Doppler shifts are significantly 
larger than the blueshifts of absorption lines from low- 
ionization-state gas. These blueshifts unambiguously 
mark outflowing gas because the absorbing gas must 
lie on the near side of the continuum source. Among 
ULIRGs, absorption in the Na I 5890, 96 doublet is 
typically blueshifted only 200-300 km s -1 at the deep¬ 
est part of the ab sorption trough dMartiul 120051 120061 
iR.unke et al.ll2005h . 

Many of our COS targets have weaker Na I absorption 
than the average ULIRG because we observed the UV 
brightest ones. Figure [12] illustrates the strong corre¬ 
lation between the Na I equivalent width and the color 
excess. We interpret this correlation as an ionization ef¬ 
fect. Sodium must be shielded from near-UV photons to 
remain ionized; it takes just 5.1 eV to remove its outer 3s 
electron. The UV selection picks the ULIRGs with less 
reddened sightlines, and sodium is therefore more likely 
to be ionized. 

The absence of Na I absorption does not indicate the 
absence of an outflow of low-ionization-state gas. The 
ionization state of sodium does not uniquely determine 
whether hydro gen is neutr al or ionized along these sight¬ 
lines dMurrav et al.l l2007ll . The wider range of ioniza¬ 
tion states accessible with far-UV transitions proves this 
point. We saw this in both Figure 0] and Figure 0 and 
we demonstrate it for the HH-dominated ULIRG with 
the best continuum S/N ratio in Figure fl3l 

The spectrum of IRAS16487+5447 also demonstrates 
how the gas kinematics changes with ionization state. 
The absorption troughs from low-ionization-state gas 
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Fig. 9.— Comparison of Lyman series and optical line profiles. The Lya absorption troughs (top row) line up with the velocities of the 
maxima in the optical emission lines. Comparison of the blended Ho; and [N II] AA6548, 6584 lines with the H (3 line shows the contribution 
of the broad line wings to the blend. The line profiles for [O III] A5007 demonstrate that broad wings are present on the forbidden lines 
as well as the recombination lines. The systemic velocities were obtained from either CO observations or the average of many optical 
emission lines as summarized in Table [T] Blended lines from the galaxy are indicated by solid, vertical lines, and Galactic absorption lines 
are indicated by dotted, vertical lines. 
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Fig. 10.— Velocity of blue wing in Lya vs. [O III] at 20% 
of the net profile area. We give the Doppler shift relative to 
the systemic velocity rather than the median velocity of the net 
emission. The plot omits the two ULIRGS without Lya detec¬ 
tions and IRAS09583+4714-B which shows a P Cygni profile (the 
blueshifted line is purely in absorption). The correlation coefficient 
is r$ = 0.93, and the null hypothesis is rejected at the 2.6a level. 
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Fig. 11. - The Lya escape fraction increases with the blueshift of 
the Lya emission wing, defined here as the velocity beyond which 
10% of the Lya emission line flux is measured. This correlation 
(r s = 0.95 and null hypothesis rejected at 2.7a) is slightly stronger 
and more significant that that of f eS c{ Lya) with the velocity of the 
Lya blue wing. 


have their deepest point near the systemic velocity. The 
greater Doppler shift of their blue wings (relative to the 
red wings) indicates a contribution from an unresolved 
outflow component, but these wings do not reach the 
-500 km s” 1 blueshift of the Lya emission wing. Ab¬ 
sorption in the medium ionization state gas probed by 
Si III probably reaches —500 km s -1 . The warmer gas 
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Fig. 12.— The COS ULIRGs \ 
cesses show the strongest Na I absorption (e.g., IRAS23365-I-3604, 
IRAS10378+1109, and IRAS0803+5243). We interpret this corre¬ 
lation as evidence for an increa sing neutral s o dium fraction wit h 
increasing dust column following lMartini (120051 ); Chen et ali i (2010 1. 
Choosing the UV-brightest ULIRGs for the COS observations se¬ 
lected the systems with the least reddening and therefore the weak¬ 
est Na I absorption. Weak Na I absorption in many of the COS 
ULIRGs does not preclude outflows of low-ionization-state gas. 


detected via O VI absorption is offset blueward of the 
systemic velocity and extends to velocities as large as 
the blueshifted Lya emission from this object. While 
general conclusions should not be drawn from a single 
object, we note that FUSE spectroscopy of nearby star- 
burst galaxies has previously revealed higher mean out¬ 
flow velocities in O VI than low-ionization-state lines in 
the same spectra dGrimes et all 12009 '). We add to the 
emerging picture the result that the Doppler shift of the 
broad emission-line wings appears to be more closely re¬ 
lated to the Doppler shift of the O VI absorption than to 
the absorption from low-ionization-state gas. 


5. ORIGIN OF BROAD EMISSION 

The presence of broad line wings on the optical emis¬ 
sion lines as well as Lya suggests gas moving at high ve¬ 
locity emits the broad, blueshifted Lya line wings. The 
predominance of blueshifted line wings over redslrifted 
wings rules out a scenario involving large random mo¬ 
tions and requires bulk flows to generate the asymmetry. 

As recognized pre vious ly in the context of the narrow 
line region of AGN (jWhittlel 11985( 1. however, the direc¬ 
tion of radial flow cannot be determined unambiguously 
from the asymmetry of an emission-line profile alone. Ei¬ 
ther wing of the intrinsic line profile can be suppressed in 
an outflow. Suppression of the red wing corresponds to 
a source of line opacity distributed throughout the out¬ 
flow region. In contrast, line opacity internal to discrete 
clouds preferentially blocks radiation from the near side 
of the outflow because the ionized side of each cloud faces 
the central source; and the resulting line profile has a red 
wing. _ _ 

We follow iSoto fe Martini (I2012T) and argue the bulk 
flow is almost certainly outward. Our primary consid- 
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Fig. 13.— Absorption lines in the COS G130M and ESI (top left only) spectra of IRAS16487+5447. Note the lack of Nal absorption 
compared to the strong Si II 1190, 1193, and 1260 lines. The Si II lines are not as broad as the Si III 1206 absorption. The O VI 
1032 absorption trough reaches even higher blueshifts of at least -500 km s 1 . (Note the blended O VI 1038 trough precludes the same 
measurement.) The O VI absorption may have some contribution from stellar winds as well as warm gas in a galactic outflow. Precise 
measurements of these troughs are sensitive to the modeling of the line spread function and the continuum placement. 


eration is the bluesliifted absorption seen in resonance 
lines such as Na I, Si II, Si III, and O VI. The gas de¬ 
tected in absorption is clearly outflowing because it must 
lie between us and the galaxy. Sources of energy and mo¬ 
mentum for the outflows include massive stars and the 
central AGN whereas, in contrast, we do not have a phys¬ 
ical reason to expect gas to be flowing into galaxies at 
speeds around 1000 km s . In the outflow interpreta¬ 
tion, emission from the fast outflow on the far side of the 
galaxy produces broad, redshifted line wings. We inter¬ 
pret the absence of such wings in our spectra as evidence 
that the receeding outflow lies largely behind the dusty 
galaxy and note that the broad, symmetric line wing s on 
the CO(1-0) profile of IRAS23365+3604 (|Cicone et al.1 

Ml . 

Our analysis generates two primary puzzles. First, 
what is the relationship between the fast outflow emit¬ 
ting the broad, emission-line wings and the lower Doppler 
shifts measured for the low-ionization-state absorption 
lines? Second, why do only the most extreme galactic 
outflows show these broad, emission-line wings? For ex¬ 
ample, the optical spectra of some high-redshift galax¬ 


ies dNewman et all l2012t iForster Schreiber et al.1 12014 
iGenzel et al .1120141 show broad, emission-line wings sim¬ 
ilar to those in ULIRG spectra. We note that one com¬ 
mon property of these high-redshift galaxies and ULIRGs 
appears to be exceptionally high SFR surface densities. 

5.1. Cool Gas Condensing Out of a Hot Wind 

We are suggesting a new picture wherein the emission 
wings come from filaments that condense out of the cool¬ 
ing hot wind. The large velocity of the blueshifted nebu¬ 
lar emission provides the primary motivation for consid¬ 
ering a power source associated with the hot wind fluid. 
Support for this picture comes from the scaling of the 
Lya luminosity with the SFR and the structure of sim¬ 
ple models for galactic winds launched from regions of 
extremely high SFR surface density. 

5.1.1. Velocity of Hot Wind Fluid 

When some of the supernova energy associated with 
star formation is ther malized, a hot wind is dri ven out of 
the starburst region. iChevalier fe Clegel (|1985l ) provided 
a simple analytic solution for the Mach number, ther- 
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mal pressure, and density in a fast, adiabatic wind. The 
wind reaches Mach 1 at the launch radius Rq where the 
mass and energy injection from the starburst region ter¬ 
minates. At r = Rq, the wind temperature in t he CC85 
solution, as shown for example in Eqn. 6 of [Zhang et al.1 

(f2f)Tl . is 

T(r) = 1.8 x 10 7 K M (£) , (5) 

where /r is the mean molecular weight per gas particle, 
and a and /3 are parameters of order unity which de¬ 
scribe the efficiency of thermalizing the supernova en¬ 
ergy (a = Ehot/E sn) and entraining gas into the wind 
(/3 = Mhot/SFR). The specific thermal energy of this 
gas is sufficient to accelerate every particle up to a ve¬ 
locity 

^=950tms-'^ 18x r i0TK . (6) 

The similarity between the specific thermal energy of 
hot winds and the specific kinetic energy of the gas emit¬ 
ting the nebular line wings motivates a conjecture that 
the highly blueshifted nebular emission arises from gas 
condensing out of a hot wind. Since the initial velocity 
of these clouds would be the speed of the hot wind, this 
picture explains their high velocities without introducing 
the problem of accelerating cool clouds to velocities much 
larger than their internal sound speeds; clouds tend to be 
destroyed by such processes. 

5.1.2. Gas Density & Instability of a Cooling Hot Wind 

iZhang et al. (|2014f ) also show that the gas density in 
the CC85 solution increases lineary with the SFR sur¬ 
face density, hereafter X* = SFR/ttRq. Adopting the 
radiative cooling fu ncti on for ga s in collisional ionization 
equilibrium (e.g.. iDraine 1201 1) and references therein), 
one finds that the cooling timescale for the hot wind, 

tcooi(R) = 7--, (7) 

(7 - 1)A (T)n e n H 

declines as the starburst becomes more concentrated, i.e., 

tcool OC S jj. 

We can define a cooling radius where the cooling time 
is reduced to the outflow timescale, t out = R/v(R), and 
this cooling radius moves inward as R coo i oc S ” 1 for the 
CC85 solution and standard cooling function. Plugging 
in typical numbers for ULIRGS, we expect the hot wind 
to be cooling very close to the starburst region, in the 
central few kpc of ULIRGs. We propose that the ULIRG 
winds start to cool at smaller radii than do winds in star- 
bursts like M82 because the higher SFR surface density 
drives winds with higher densities and smaller cooling 
radii. 

The high gas density at relatively small radii leads 
to a much higher emission measure at the cooling ra¬ 
dius than would be expected in more typical starburst 
galaxies. As the wind cools, we expect a highly unsta¬ 
ble situation to develop as the temperature drops from 
the free-free portion of the cooling curve to a minimum 
(near T ss 2 x 10' K for solar metallicity gas in CIE). 
Between 10 5 < T(K ) < 10 73 , the cooling rate per atom 


begins t o increase rapidly as the gas cools, A(T) oc T 0 7 
(|Drainel |20Tll). Resolving this process and the associated 
condensations of clouds in numerical simulations would 
be enlightening. 

5.1.3. The Correlation of Lya and Bolometric Luminosity 

The ULIRG Lya luminosities span a much larger dy¬ 
namic range, about four orders of magnitude, than do 
their far-infrared luminosities. Figure IT4l illustrates that 
a power-law scaling, L(Lya) oc L 2 FIR , describes the lumi¬ 
nosities of the ULIRGs not dominated by AGN. Includ¬ 
ing the two AGN, IRAS11598-0112 and IRAS01003-2238, 
the exponent steepens to approximately four. 
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Fig. 14. — The ULIRGs show a much larger dynamic range in Lya 
luminosity than they span in far-infrared luminosity. Excluding the 
two non-detections (triangles), we find a strong, positive correlation 
between the Lya luminosity and the far-infrared luminosity. The 
dashed lines illustrate a power-law scaling, L( Lya) oc L 2 R . 


A quadratic scaling like this has been predicted, but 
not observed, for the X-ray luminosity emitted by star- 
burst winds. Using the analytic solution for a strong 
wind driven out of a starburst region w here most th e su¬ 
per no va energy has been t hermalized (jGhevalier &; Cleed 
I1985T) . IZhang et al.1 (1201411 calculate the 0.5-8.0 keV lu¬ 
minosity of the hot wind fluid and find Lx oc SFR 2 . 
This scaling requires constant thermalization and mass- 
loading efficiencies. It is interesting that IZhang et al.1 
(|2014ll found that the Lx — SFR relation measured 
among star-forming galaxies is shallower than their pre¬ 
diction. The shallow relation can be attributed at least 
in part to the diffuse X-ray emission from galaxies being 
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brighter than the hot wind fluid. The hot wind fluid is 
extremely challenging to detect directly due to its low 
emission measure. The extended X-ray emission sur¬ 
rounding starburst galaxies is widely believed to emanate 
from intermediate temperature gas at the boundary be¬ 
tween gas ph ases because it is spati ally correlated with 
Ha emission dStrickland et al.l 12004) . That the Lya lu¬ 
minosity increases as roughly the square of the SFR, as 
traced by Lfir , suggests to us that this luminosity may 
b e physic ally rel ated to the cooling of the hot wind. 

iZhang et al.l (|2014l ) predict the 0.5-8.0 keV X-ray lu¬ 
minosity of the hot wind fluid per unit SFR, 


L^ ot [0.5 


8.0 keV] = 10 43 ergs s" 1 



where Rq descr ibes the radius of the starburst region. We 
substitute the iKennicuttl (|1998T ) relation between SFR 
and far-infrared luminosity and obtain 


T hot 
L X 

Lfir 


0.0045 


/ SFR \ 
\100 M 0 yr-V 


/ 200 pc\ 

V Ro )' 


(9) 


The dashed lines in Figure Q3] compare the Lya lumi¬ 
nosity to L the predicted cooling luminosity. Among 
the starburst-dominated ULIRGs, the Lya luminosity is 
~ 0.3% of the pedicted cooling radiation from the hot 
wind fluid. The closet analogy may be the cooling of hot 
gas in the centers of galaxy clusters where calculations 
for T ~ 10 7 K gas cooling non-uniformly show of order 
1% o f the losses in strong UV and optical emission lines 
(jVoit et al.l 1 19941b 

As shown in Sect ion l4~Tl most of the Lya photons emit¬ 
ted near the systemic redshift are scattered by neutral 
hydrogen; but the blueshifted wing on the Lya profiles 
does not arise from scattering and may therefore be emit¬ 
ted at larger radii. The luminosity of just this highly 
blueshifted component may be more directly related to 
the cooling wind than is the total Lya luminosity. 

We compare just the Lya luminosity emitted blue- 
ward of —500 km s -1 with the far-infrared luminosity 
in Figure 1151 The luminosity of the highly blueshifted 
Lya emission increases with the far-infrared luminos¬ 
ity. Among the HH-dominated ULIRGs, the highly 
blueshifted Lya emission is about 0.1% of the far-infared 
luminosity. The blueshifted Lya emission is a larger 
fraction of Lfir , roughly 5%, in the AGN-dominated 
ULIRGs. The steeper than quadratic rise in y a 
with Lfir between the HH-dominated AGN-dominated 
ULIRGs shows that the emergence of the AGN increases 
the Lya emission from ULIRGs. 


5.2. Fermi Acceleration of Lya Photons 

The framework for the fast outflows should explain 
the shock speeds measured from optical emission-line 
ratios in addition to the lower Doppler shifts of res¬ 
onance absorption lines. The scattering of Lya pho- 
tons across shock fronts could be one such mechanism. 
iNeufeld fe McKee] (|1988l) have shown that scattering of 
resonance photons across shock fronts systematically 
blueshifts the photons. Since the photons pick up a 
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Fig. 15.— The Lya luminosity in the blue wing, defined here as 
v < —500 km s“ 1 increases with the far-infrared lu mino sity in a 
manner similar to the total Lya luminosity in Figure [Til 

Doppler shift of 0.51^ on each scattering, the blueshift 
exceeds the shock velocity after only a few reflections 
across the shock front. 

We considered the shocks identified by iSoto et all 
(|2012f) IRAS00262+4251 as a preliminary test of this 
idea. A shock speed of v s h = 225 ± 50 km s _1 fit 
the optical emission-line ratios well. In the Fermi ac¬ 
celeration model, only three such scattering encoun¬ 
ters are sufficient to produce the observed line wing in 
IRAS00262+4251 which has Vs 20 = -326 km s" 1 . The 
Doppler shift of the peak scattered emission scales as 
= 270(iV 2 ou s7 ) 1 / 3 in their model. 

The scaling of this model predicts the velocities of the 
blue wings increase with shock speed and the column 
density of the neutral gas, N 20 , Comparison of the shock 
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velocities from iSoto fe Martini (l2012f) . reproduced here 
for convenience in Table 0 against the velocity of the 
blue wings from Table [4] reveals no correlation between 
the two velocities. It follow that the column density of 
neutral gas would have to be correlated with the velocity 
of the blue wing for the data to support Fermi accel¬ 
eration acceleration as the primary explanation for the 
Lya line wings. We have argued, however, that the col¬ 
umn densities are lowest in the objects with the bluest 
line wings because we measure the largest f esc (Lya) in 
those objects. Hence, we find no compelling evidence 
that Fermi acceleration shapes the blue wings on the 
ULIRG Lya profiles. We emphasize, however, that the 
correspondence between the shapes of the Lya and op¬ 
tical line wings varies considerably among the ULIRGs. 
And Fermi acceleration seems like an important physical 
process that could play a role in some objects although 
it awaits compelling empirical support. 


5.3. Relation between Outflow Velocity and Galaxy Size 

Some wind models predict the fastest outfl ows from the 
most compact starbursts. For example, iMurrav et al.1 
<f2on argue that the cluster luminosity cannot grow 
much beyond its Eddington limited value due to outflows 
driven by radiation-pressure which disrupt giant molecu¬ 
lar clouds faster than the cloud dynamical timescale. Mo¬ 
mentum conservation requires terminal velocities scale 
as 

Vterm OC Rq 1/2 , (10) 


where Rq again describes the radius of the region where 
the wind is launched. In general support of this scenario, 
the galaxies with the most blueslrifted resonance absorp¬ 
tion turn out to indeed be very comp a ct with radii of 
rough l y 100 pc dTremonti et al.1 (120071) : iHeckman et al.1 
(|2011h : iDiamond-Stanic et al.1 ( 2012D ). The discovery of 
fast outflows from ULIRGs therefore begs the question 
of the starburst size in ULIRGs and the SFR surface 
density. 

The starbursts in ULIRGs have long been considered 
concentrated; and one of the main differences between 
ULIRGs and the sub-mm-selected galaxies at high red- 
shift is the smaller size of the starburst region in the 
former. One estimates the effective size of the starburst 
region in a ULIRG from a black body argument. This 
approach assumes the galaxy is optically thick. For ex¬ 
ample, iScovillel (l2013f) give 


Rfir = 1-09 kpc 


Lfir 
10 12 L 0 


1/2 


(¥)'■ 


( 11 ) 


where the fiducial dust temperature, Td, was chosen to 
represent starburst dominated ULIRGs. If the source is 
optically thin, then this equation provides a lower limit 
on the radius. 

Substituting the far-infrared luminosities from Table |Tj 
into Eqn. El we estimate the starburst radii range from 
0.85 to 1.2 kpc and greatly exceed those of the DCO 
and the massive, post-starburst samples. This result is 
clearly sensitive, however, to the assumed dust temper¬ 
ature. If Td ~ 100 K in the AGN-dominated ULIRGs 
IRAS11598-0112, then Eqn. fill gives a size for the emit¬ 
ting region of just Rfir » 0.18 kpc (100 K /Td) 2 . 


From Eqn. El we can see that far-IR sizes for 
starburst-dominated ULIRGs will increase as the square- 
root of their luminosity, Rfir oc L rir . Hence, any cor¬ 
relations between starburst size and fast outflows in the 
COS ULIRGs sample reduces to the strength of the corre¬ 
lation between the blue wings and the far-infrared lumi¬ 
nosity. Regardless of whether we consider Lya or optical 
emission lines, the line wings extend to higher velocities 
in more luminous systems 03 We conclude that among 
the starburst ULIRGs, those with larger far-IR sizes tend 
to show more extended blue wings in Lya. Only when 
the starburst dominated ULIRGs are compared to the 
AGN dominated ULIRGs as a class would we recover 
the inverse scaling of faster outflows from smaller launch 
regions that the escape velocity argument would favor. 

For completeness, we also estimate NUV sizes using 
the WFC3/F225W images shown in the middle panel 
of Figure [3] Specifically, we measured the area of the 
aperture enclosing all the pixels with surface bright¬ 
ness greater than a fiducial value of E(/ llt (A2366) = 
7.1 x 10 26 ergs s -1 cm -2 Hz -1 kpc -2 . This threshold 
represents a SFR of 0.1 M 0 yr -1 kpc -2 for continuous 
star formation (of duration greater than 100 Myr with 
a Salpeter initial mass function) and is of interest be¬ 
cause it is widely believed to be the minimum value re¬ 
quired to drive galactic winds IHeckman] (|2003l ). We find 
values of this Rnuv = (A/tt) 1 / 2 from 0.5 to 8.2 kpc. 
We find no correlation between these NUV sizes and 
the FIR sizes; the assumed dust temperature is irrele¬ 
vant to this comparison because our imaging sample ex¬ 
cludes the ULIRGs with the strongest AGN. We find that 
IRAS12071-0444 has a large size Rnuv = 5.3 kpc which 
supports our argument lSection l2.1.2[) that the starburst 
is at least as luminous as the AGN in this object. These 
NUV sizes show no correlation with Vb 20 - 

The notable point here is that for a constant dust tem- 
perture, the larger starbursts show higher outflow veloc¬ 
ities, and this trend is opposite to both the theoretical 
prediction for clouds accelerated by radiation pressure 
and the Doppler shifts of resonance absorption lines in 
spectra of extremely compact starbursts. In our pro¬ 
posed picture, however, the broad wings on emission lines 
come from condensation of the hot wind whereas the 
resonance absorption lines come from clouds accelerated 
by some combination of radiation pressure, cosmic rays, 
and/or the hot wind. The difference between the scal¬ 
ing relation for the broad line wings and the resonance 
absorption lines therefore poses no obvious contradiction 
but should be explored further via numerical simulations. 
We also emphasize the need for resolved measurements of 
the starburst sizes in ULIRGs and/or SED constrained 
dust temperatures; short of these, the size estimates too 

11 Formally, the correlation coefficient between Vg 20 ,sys (Lya) 
and L(Lya) is Rg = —0.62 at only 1.74(7 significance. The Lya lu¬ 
minosity correlates more strongly with the velocity of the blue wing 
on the [O III] line, VB20,sys ([0//J]); and we measure a correlation 
Rg = —0.80 at 2.26cr significance. The correlation of the Lya lumi¬ 
nosities with the far-infared luminosities is not one-to-one; and it 
is interesting that the correlations weaken when we replace L(Lya) 
with Lfir • The correlation coefficient between VA 20 ,si/s(Lya) and 
Lrir ' s Rs = —0.60 at 1.79a significance. The far-infared lumi¬ 
nosity correlates less strongly with the velocity of the blue wing 
on the [O III] line, Vr 20 ,sys([OIII]); and we measure a correlation 
Rg = —0.57 at 1.61(7 significance. 
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uncertain to consider them compelling evidence for the 
condensation scenario. 

5.4. Evolution of Lya with Merger Progression 

ULIRGs are examples of L* galaxies caught during 
their transition to the red sequence. As we do not yet 
understand what makes, and keeps, galaxies ’red and 
dead,’ direct measurements of feedback in these objects 
are of particular interest. What we know already is that 
towards the final merger stage a larger fraction of the 
bolometric luminosity in ULIRGs appears to be powered 
by an accreting, supermassive black hole (lYeilleux et al.1 
1200911 . and the accompanying increase in bolometric lu¬ 
minosity correlates with an increase in the speed of the 
molecular outflows (lYeilleux et alJ2013j l. It is interesting 
to explore how the broad line wings fit into this picture. 

To make this comparison, we refer back to the 
AGN fractions introduced in Section 12.1.21 and supple¬ 
ment those classifications with measurements of the far- 
infrared luminosities. In Figure fill for example, we saw 
that the ULIRG with the dominant AGN is 25% more lu¬ 
minous than any of the other objects in the sample; this 
result is consistent with the trend that more luminous 
ULIRGs have a higher AGN fraction on average. What 
Figure [TT] shows that had not been seen before is a high 
ratio of Lya to far-infrared luminosity in IRAS 11598- 
0112 and in the ULIRG with the next highest AGN frac¬ 
tion, IRAS01003-2238. We note that the enormous in¬ 
crease in the strength of the Lya emission (relative to 
Lfir in these two ULIRGs) and attribute it to the emer¬ 
gence of the AGN. This increase in Lya production may 
be directly related to the AGN spectrum or (indirectly) 
to changes in the ISM induced by AGN feedback which 
enhance the escape fraction of Lya photons. The take¬ 
away point is that the ULIRGs with a dominant AGN 
show evidence for additional production of Lya emission. 

Starbursts, in contrast to AGN, produce the majority 
of the bolometric luminosity in the typical COS ULIRG, 
so an important question is whether the blueshifted wing 
on the Lya profile in a starburst-dominated spectrum 
could be directly related to an AGN, albeit a less lu¬ 
minous one. In our opinion, the properties of the star- 
burst ULIRGs measured to date suggest the answer is no. 
First, we have demonstrated that the kinematics of the 
Lya wings indicate they come from the same parcels of 
gas emitting the weaker, but equally broad, wings on the 
optical emission lines. Among the starburst-dominated 
ULIRGs, the region of broad optical emission subtends 
an angle along the ESI slit that corresponds to up to 
20 kpc with an average radius of 3 kpc. Due to the 
large size of the excited region, it is more natural to at¬ 
tribute the shock-like line ratios of the broad component 
to out flows rather than to the narrow-line re gion of an 
AGN dSoto et al.ll2012t iSoto fe Martini I2012H . Second, 
we found that the ratio L(Lya)/ Lrir is over an order of 
magnitude higher in the AGN-dominated ULIRGs than 
in the starburst-dominated ULIRGs. This abrupt jump 
indicates strong evolution in either the production mech¬ 
anism or escape fraction of Lya emission as the AGN 
luminosity grows. We would have expected a contin¬ 
uous increase in L(Lya)/L rir with Lrir if the AGN 
produced the Lya emission in the starburst-dominated 
ULIRGs. 


In our initial discussion of Figure [3] we noted that the 
Lya wings extend to larger bluesliifts as the Lya lumi¬ 
nosity increases. Since the AGN fraction is known to 
grow with increasing bolometric luminosity, we plot the 
velocity of the blue wing against far-infrared luminosity 
in Figure 1161 In this plane, the correlation between the 
Doppler shift and luminosity is less pronounced (than 
Vb 20 vs. L(Lya)) because IRAS01003-2238 has the sec¬ 
ond highest Lya luminosity but an average Lrir ■ It is 
the broad Lya and optical emission-line wings on the 
mini-BAL, IRAS11598-0112, which give the impression 
of increasing Doppler shift with increasing Lrir in Fig- 
u re m 

lYeilleux et al.l (|2013ll found that the velocities of the 
molecular outflows in ULIRGs increased gradually with 
AGN fraction but showed no significant dependence on 
the SFR. While the overall correlation seen in FigurefTfilis 
not inconsistent with the trend seen among the molecular 
outflows, neither do the data strongly support it when 
one looks closely at the AGN fraction among individual 
systems. For example, the Seyfert 2 galaxies IRAS01003- 
2238 and IRAS12071-0444 have AGN which produce up 
to 50% of LroI , but their Vb 20 measurement are nearly 
the same as that of IRAS10378+1109 and lower than 
that of IRAS23365+3604, two ULIRGs with small AGN 
fractions. 

Detailed mapping of the Lya, optical, and molecular 
line wings is needed to clarify the relationships among 
these outflow dia gnostics. As a prelude to such work, we 
simply note that lYeilleux et all (I2013f) provide measure¬ 
ments of the molecular outflow speeds for two of the COS 
ULIRGs: IRAS12071-0444, V M {OH) = -321 km s" 1 , 
and IRAS23365+3604, U 84 (OH) = -604 km s" 1 . The 
blueshifts of the broad Lya wings (VB 2 o(Lya) = - 
743 km s -1 and -1046 km s -1 , respectively) indicate 
faster outflows than do the OH measurements. Whether 
these discrepancies reveal different components of the 
bulk flow or merely variations in measurement technique 
remain to be determined. 

Further exploration of the variation in outflow diag¬ 
nostics along the merger sequence promises to provide an 
empirical description of how feedback disrupts the ISM 
of a galaxy. We will show in Section [6] that modeling the 
Lya wings places significant limitations on the column 
density of neutral hydrogen along the ULIRG sightlines. 
In the most extreme cases, which are the AGN, the feed¬ 
back has generated holes in the ISM which have become 
optically thin to Lyman continuum radiation. 

6. RADIATIVE TRANSFER MODELING 

To better determine the physical conditions along the 
UV-selected sightlines, we performed Monte-Carlo ra¬ 
diative transfer simulations for the Lya emission. We 
first tried a simple model consisting of &_ thin sp heri¬ 
cal shell moving radial ly at velocity V s (lAhn et al.ll2003fc 
lYerhamme et alj|20 06). The shell has a column density 
N(HI) in neutral hydrogen and dust optical depth 
We found it difficult to produce emission at the systemic 
velocity and the broad blue wings with this model, so we 
added a single parameter, fc, to the shell model to de¬ 
scribe the covering factor of neutral gas; in this fiducial 
model, a fraction of sightlines 1 — fc = flhoie/47r in- 
tersect a ‘hole ’ whic h contains no neutral gas and dust 
S Behrens et al\\20lA) . 
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Fig. 16.— Doppler shift of Lya blue wing increases with ULIRG 
luminosity. The velocity measurement for IRAS 11598-0112 is an 
upper limit because the wings of the Lya: profile are blended with 
metal lines; comparison to the unblended wings on the optical emis¬ 
sion lines suggests, nonetheless, that the wings on the deblended 
Lya: are the broadest in the COS ULIRGs sample. When the 
ULIRGs with dominant AGN (IRAS11598-0112, IRAS01003-2238, 
and possibly IRAS 12071-0444) are removed, the tendency for the 
velocities to increase with Lrir is much less pronounced. 

The spectrum of Lya photons emerging from a spher- 
ical shell that contains a hole depends on direction (see 
Behrens et al. 2014 for a detailed analysis). Here, we 
capture this effect in a simplified (but faster) way, which 
is based on the fact that sightlines through the holes con¬ 
tain a superposition of the intrinsic (unscattered) and the 
scattered spectrum. We denote the (normalized) spec¬ 
trum as a function of velocity offset v (of Lya photons 
that scatter at least once with 5 sca t{v)) while we denote 
the intrinsic spectrum with S'int(h). Under the assump¬ 
tion that scattered radiation emerges isotropicalhQ the 
normalized spectrum of Lya photons along sight lines 
passing through holes is 

b'holes U) OC /"scat $scat (c) “t“ *Sint (c) , (12) 

in which / scat denotes the fraction of photons that scat¬ 
ter, which depends on the hole covering factor 1 — fc- 

We carried out th e calculations using t he code de¬ 
scribed previously in iDiikstra et al.1 (2006a.b). We used 
10 5 Lya photons initially distributed in frequency ac¬ 
cording to a model line profile matched to the H/3 pro¬ 
file. This intrinsic Lya spectrum implicitly associates 
the source of the Lya photons with recombination in 
H II regions. We start with a Doppler parameter 
b = 12.9 km s -1 but find it necessary to increase it to 
40-50 km s -1 to fit IRAS01003-2 238, IRAS10378+1109 
and IRAS09584+4714-B Since iBehrens et al.1 (|2014T l 
showed that the spectrum which emerges from a shell 

12 Under this assumption the scattered component can be viewed 
as a single isotropically emitting source with a total flux /scat Ftot 
(where Ftot denotes the total emitted flux). We only have an unob¬ 
scured view of the central source along sight lines passing through 
holes. This unobscured source simply appears as an isotropic 
source with a total flux Ftot- Therefore, the total flux along sight 
lines with holes is (/ sc at + l)-Ftot- 


containing holes depends on direction, we also explored 
the directional dependence of the spectrum. We studied 
Lya transfer through shells that contain holes and pre¬ 
dicted spectra for those sightlines passing through the 
holes. We found the best fitting models were consistent 
(within a factor of roughly two) with those obtained from 
the fiducial model. We use our fiducial model, however, 
because it requires significantly fewer photons to con¬ 
struct spectra (it assumes spherical symmetry), and is 
therefore better suited to exploring the parameter space 
describing the shell. 

Our approach was to find shell models that reproduce 
the basic properties of the COS Lya profiles and then 
gain physical insight from the model parameters. The 
detailed Lya shape will clearly depend on factors in ad¬ 
dition to the model parameters including viewing angle, 
source size, and shell thickness which are not directly 
constrained. Figure [17] shows examples of models that 
describe the main features of each of the 9 COS Lya 
profiles. The parameters, which are shown on each plot, 
span a substantial range but generally indicate low col¬ 
umn densities of neutral gas (17.0 < log Nhi{ cm -2 ) < 
20.3), median shell velocity of 80 km s -1 , covering frac¬ 
tions from 0.3 to 1, and relatively low dust optical depths 
{jd < 2). The model parameters, however, are likely rep¬ 
resentative of the true values (along the least obscured 
sightlines) to within a factor of a few as supported, for 
example, by our finding similar parameters for models 
in which we accounted for the directional dependence of 
the spectra that is naturally associated with shells that 
contain holes. 

The absorption troughs in the COS spectra deter¬ 
mine the velocity of the neutral gas in the shell model. 
In our description of IRAS09583+4714-B, for example, 
an expanding shell absorbs the Lya photons emitted 
slightly blueward of line center due to the Doppler shift 
of the source in the frame of the shell. Figure [171 shows 
the fit with log W//(cm -2 ) = 20.0, V s = — 80 km s -1 
(blueshift), and near unity covering fraction, fc = 0.9. 
The flux minimum in the IRAS23365+3604 Lya pro¬ 
file is also well modeled with a shell moving, but the 
strong transmission at the systemic velocity precludes 
neutral gas at the systemic velocity and requires a shell 
velocity near 400 km s -1 . Similarly, a shell moving 
at±150 km s -1 allows transmission near the systemic 
velocity in IRAS11598-0112. The redshifted absorp¬ 
tion in the IRAS16474+3430 COS spectrum, which ex¬ 
hibits an inverted P Cygni profile with blueshifted emis¬ 
sion, is well described by an shell falling towards the 
galaxy at roughly 250 km s -1 . We produced the red- 
shifted dip in the IRAS01003-2238 Lya profile with in¬ 
fall at 200 km s -1 . In contrast to these neutral gas 
flows, a static shell or simply the galactic ISM absorbs 
Lya photons near the systemic velocity in the spec¬ 
tra of IRAS00262+4151 and IRAS12071-0444; and the 
low shell velocities in IRAS10378+1109 (30 km s -1 ) 
and IRAS 16487+5447 (50 km s -1 ) could plausibly arise 
from wavelength calibration errors rather than bulk mo¬ 
tion. 

Transmission at the shell velocity can be recognized 
by emission lines within saturated absorption troughs or 
absorption troughs that are not black and requires ei¬ 
ther partial covering, and therefore some leakage of the 
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Fig. 17.— Radiative transfer models for Lya profiles. Starting from an intrinsic line profile (dashed line), we model the transfer of Lya 
photons through a shell to obtain the emergent spectrum (red curve). These models provide a rough description of the observed Lya line 
profiles (blue line), and Table [6] lists the parameters of the model shown for each target. The sign convention used here for purposes of 
display defines redshifts as positive velocities. 


intrinsic spectr um, or a low column N(HI) < 10 17 cm~ 2 
(jBehrens et all[20 1 11 ). The model for the Lya spectrum 
of IRAS01003-2238 requires columns of neutral hydro¬ 
gen which are optically thin to the Lyman continuum. 
The transmision at zero velocity in the IRAS11598-0112 
COS spectrum required a model with a low H I col¬ 
umn to N(HI ) < 10 18 ' 8 cm -2 . It is interesting that 
while IRAS23365+3604 has substantial emission at the 
systemic redshift, Lyman continuum radiation will be 
absorbed by the shell at 400 km s" 1 and will not es¬ 
cape the galaxy. These results are consistent with a 
decline in the neutral gas column as the merger pro¬ 
gresses since IRAS23365+3604, IRAS11598-0112, and 
IRAS01003-2238 are classified, respectively, as a close 
binary/pre-merger, compact merger, and old merger. 

The dust optical depth in the models can be directly 
compared to the measured nebular reddening. The dust 
optical dept h near Lya is = 10.0 E(B — V) for the 
iCardelli et all (11989fl reddening curve, so it is trivial to 
convert the color excess listed in Table [5] to a meaured 
dust optical depth. Comparison to the model in Ta¬ 
ble [6] shows that the dust optical depth inferred along 
the UV sightlines is generally lower than that measured 
towards the HII regions. A number of physical interpre¬ 


tations are possible, and we note only that this result 
suggests the UV and optical sightlines do not sample ex¬ 
actly the same regions of these complex systems. 

This modeling exercise also provides insight into the 
origin of the broad wings of the COS Lya profiles. The 
blueshifted emission in the Lya line can be attributed 
primarily to the blue wing of the intrinsic profile. 

6.1. Comparison to Previous Modeling 

lYerhamme et al.l (|2014fl and lBehrens et ahl (|2014f l have 
also recently explored the Lya profiles of shell models 
with low neutral columns and partial covering. Our re¬ 
sults strongly support the primary conclusion of those 
studies, namely that the emergence of Lya emission near 
the systemic velocity boosts the probability that ionizing 
photons may leak through. 

The major difference is that we start with intrinsic 
profiles which have blue wings. The intrinsic Lya profile 
may have a red wing in addition to the blue wing, but we 
start with the simplest profile with which we can model 
the data. 

The Lya signatures that betray the escape of Lyman 
continuum photons may, ho wever , depen d quite a lot on 
the model. For example. lYerhamme et al.l (1201411 suggest 
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that the peak Lya emissi on from density-bounded HII 
regions (|Osterbrock &; Ferlandl 12006!) shifts by less than 
150 km s -1 from the systemic velocity. Close inspection 
of Figure [T7] shows that the shift between the absorp¬ 
tion feature and the maximum emission ranges from 150 
to 300 km s -1 in IRAS00262+4251, IRAS10378+1109, 
IRAS16474+3430, and IRAS12071-0444; our models for 
these profiles all have logAbj/(cm -2 > 18.0), so our 
results are consistent with the hypothesis that no Ly¬ 
man continuum escapes. Our model for IRAS01003- 
2238, however, requires a low enough column den¬ 
sity that we expect a significant Lyman continuum es¬ 
cape fraction, yet the emission peaks are still shifted 
by more than 150 km s^ 1 , specifically +200 km s^ 1 
and —300 km s _1 from the absorption feature visible at 
200 km s -1 in Figure |T7] Adding an anisotropic gas 
distribution to the fiducial model could produce veloc¬ 
ity shifts as large as we measure (even for low neu- 
t ral gas colu mns) when viewed from specific directions 
(IBehrens et al.ll2014fl. 

The lVerhamme et al.1 (12014H study concludes that peak 
separation will be less than 300 km s -1 when Lyman 
continuum escapes. The peak separation in our radiative 
transfer simulation of IRAS01003-2238 is larger, close to 
500 km s -1 . In our model, we must have the shell mov¬ 
ing at 200 krn/s to get the absorption right. We cannot 
increase the neutral hydrogen column density of the shell 
much above 10 1 ' cm" 2 , or we run into trouble with the 
observed damping wing. The photons that are absorbed 
in this shell and then backscatter pick up a redshift of 1 
times the outflow velocity (because they appear practi¬ 
cally at resonance in the frame of the gas). In this model, 
the peak-tro ugh separation is twice t he shell velocity. 

Figure 3 of lVerhamme et ahl (120141 1 indicates that the 
blue peaks go away when the shell velocity exceeds 
100 km s -1 . Our model for IRAS01003-2238 (the thin 
case) as well as those for IRAS16474+3430, IRAS11598- 
0112, and IRAS23365+3604 all have blue peaks even 
though the shell velocity in the model is larger than 
100 km s -1 . This difference further emphasizes our re¬ 
sult that the intrinsic Lya profiles of the COS ULIRGs, 
prior to processing by radiative transfer, must have high 
velocity blue wings and suggests the above quantitative 
result might change with assumptions about the shape 
of the intrinsic line profile. 

6.2. Other Evidence for the Escape of Ionizing 
Radiation 

While radiative transfer models clearly provide ex¬ 
tremely useful insight into the physical conditions in 
ULIRGs which shape their unusual emergent Lya pro¬ 
files, the simple geometry inherent to these models proba¬ 
bly precludes detailed quantitative assessment of Lyman 
continuum escape from the separation of profile peaks 
alone (see § 16. .71 for a discussion). What they confirm, 
however, is that the shapes of the Lya profiles is the COS 
spectra make ULIRGs leak some Lyman continuum ra¬ 
diation. A quantitative measurement of this leakage re¬ 
quires direct detection of the Lyman continuum. Short of 
that, however, we can only consider other indirect tracers 
of the optical depth. 

In particular, the ESI spectra cover the AA3726, 29 dou¬ 
blet as well as the [O III] AA4959,5007 doublet. Ex¬ 


treme flux ratios, F[OIII\/F[OII\> 9 can be an indi¬ 
cation of density boun ded HII reg ions which have un¬ 
dersized 0 + zones dJaskot fe Oevl f2Q13h . Hence, our 
radiative transfer modeling results predict a very high 
flux ratio for IRAS01003-2238. In support of this pre¬ 
diction, we confirm that this AGN has the highest flux 
ratio in the COS ULIRG sample, F[OIII]/F[OII] « 
6.0. The only other COS ULIRG with a flux ratio 
greater than 2 is IRAS12071-0444, a possible AGN, with 
F[OIII]/F[OII ] ~ 4.5 and a much higher fitted neutral 
gas column in Table [6] These line ratios would hardly 
be compelling evidence for Lyman continuum escape on 
their own because the high ratios could also be produced 
by large ionization parameters in these nebulae, but they 
support our conclusions that (1) IRAS01003-2238 almost 
certainly has the lowest neutral gas column in the COS 
ULIRGs sample, and (2) it is plausibly optically thin to 
the Lyman continuum. 

Another approach would be to measure the relative 
residual intensity in the cores of the strongest low- 
ionization-state abs o rption lines dHeckman et al.1120011: 
Grimes et al.1 120091 IHeckman et al.1 120111 I Jones et all 
2013T) . Partial covering of low ionization state gas can be 
measured from the dept h of saturated absorption lines 
(|Martin fe Bouchl 120091) . Deeper COS spectra, such as 
one might obtain for purposes of directly measuring the 
Lyman continuum, could be used to measure the cover¬ 
ing fraction via lines of Si II or C II. 

6.3. Other Explanations for Lya Escape at Line Centre 

Shell models are frequently used to model Lya lines 
because of the relative ly small number of parameters 
(iVerhamme et all 120081 ). We found that ordinary shell 
models could not fit these lines. We therefore intro¬ 
duced a hole, which represents a simple way to include 
anistropic escape. These models allow for more Lya pho¬ 
tons to escape without scattering, and at line centre. We 
stress that however, that the escape of Lya photons at 
line center is not a unique signature of clear sight lines to 
the nebular regions. Two alternative explanations exist 
in the literature: 

1. iGaravito-Camargo et al.1 (120141) have shown that 
Lya transfer through rotating optically thick 
clouds can give rise to Lya photons escaping at 
line center. 

2. The ‘Neufeld mechanism’ dNeufeldl I1991D allows 
photons emitted at line center to escape at line 
center from an optically thick medium if they only 
scatter off the surfaces of chimps in a multiphase 
medium. In order for the Neufeld mechanism to 
transmit photons at line center, we need station¬ 
ary clumps (line center photons would not ap¬ 
pear at line center in the frame of out-or-inflowing 
clumps). We can replace our ‘empty’ holes with 
holes that are fully covered by stationary clumps 
with Ahi ~ 10 17 cm -2 . This would not affect 
the RT much in models in which shells contain 
Ahi 10 1 ' cm -2 , but it would modify the ra¬ 
diative transfer in e.g. object IRAS01003-2238 in 
which we inferred Ajji = 10 1 ' cm -2 . 

Our models for the observed Lya lines are therefore 
likely not unique: for example, it remains possible that 
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some observed line profiles allow our ‘holes’ to contain 
enough HI to be opaque to ionising radiation. However, 
in some models (e.g. IRAS01003-2238) there is so little 
scattering at all frequencies that this is likely not possi¬ 
ble. Addressing both the effects of rotation and clumps 
in a systematic way is well beyond the scope of this pa¬ 
per, but we expect to explore this in more detail in a 
separate paper, in which we will also study the addi¬ 
tional constraining power of observing the Ly/3 regime 
(repeated surface scattering of Ly/3 photons is limited by 
the conversion probability of Ly/3 into Ha + 2q contin¬ 
uum) . 

7. SUMMARY AND IMPLICATIONS 

Feedback from star-formation on interstellar gas is cen¬ 
tral to theories for galaxy formation and evolution, and 
the most plausible mechanism for powering the removal 
of gas and metals from galaxies remains a galactic wind. 
The challenge at present revolves around translating the 
observable signatures of these winds into physical prop¬ 
erties. One promising path forward would be to connect 
measurements of outflows to simple, analytic descriptions 
of the wind physics. Our confidence in this mapping 
would be enhanced if we could understand different ob¬ 
servational diagnostics in the context of a simple physical 
framework. 

In this paper, we used new far-UV spectroscopy from 
HST/COS in concert with high-resolution optical spec¬ 
troscopy to compare signatures of galactic outflows. The 
sample includes the most luminous and dusty starbursts 
found at low redshift, yet we detected Lya emission from 
9 of the 11 galaxies. Even more surprising, however, 
only the spectrum of the galaxy that is not a ULIRG 
(IRAS09583+4714-B) shows the classic P Cygni profile 
produced by resonance scattering off an expanding, neu¬ 
tral medium. We do identify signatures of resonant scat¬ 
tering in the Lya profiles; these include absorption near 
the systemic velocity, broad line cores compared to non- 
resonant lines, and a redshifted component absent in the 
optical lines. 

The unanticipated feature of the Lya line profiles is 
their broad, blue wings. The blueshift of this emission 
exceeds 500 km s -1 in five of the COS ULIRGs. It seems 
plausible, perhaps likely, that these wings arise from 
AGN activity in IRAS11598-0112 and IRAS01003-2238, 
the two objects where an active nucleus provides > 50% 
of the bolometric luminosity. In contrast, the total Lya 
luminosity is over an order of magnitude lower (when 
measured relative to the far-infrared luminosity) in the 
ULIRGs with starburst-dominated spectra and cool far- 
infrared colors. We find it more natural to consider ori¬ 
gins for the broad wings related to the starburst winds 
in these systems. 

We attribute the large Doppler shifts (of order 
1000 km s" 1 ) to bulk outflow based on (1) the close cor¬ 
respondence between the shapes of the blue wing on the 
Lya and optical emission-line profiles, (2) the similar at¬ 
tenuation of the Lya line and the underlying continuum, 
and (3) the need to include these wings on the intrinsic 
Lya profile in radiative transfer modeling of the emergent 
line profile. This blueshifted Lya emission must come 
from the near side of the galaxy and hence outflowing; 
otherwise the resonance photons would be scattered by 
the neutral hydrogen column. The velocity of the out¬ 


flowing gas in the line wings is two to three times the 
escape velocity from the starburst. 

The high velocity of the Lya emission compared to 
the lower blueshift of resonance absorption lines poses a 
paradox. The low-ionization-state gas traced by Na I and 
Si II absorption is moving outwards more slowly than the 
denser gas seen in emission. The Doppler shifts of the 
emission-line wings also exceed the shock velocities fit to 
emission-line ratios by a factor of two to three. We do 
detect gas in absorption at these high velocities via the 
O VI 1032, 1038 doublet. This tantalizing evidence for 
an association with a higher-ionization-state component 
of the outflow is currently limited to a few ULIRGs, how¬ 
ever, because the GO 12533 observations were designed 
to study Lya emission rather than UV continuum. 

Motivated by our initial goal of tying attributes of the 
line profiles to simple, physical models for galactic out¬ 
flows, we advance the conjecture that this cool, fast gas is 
a direct prod uct of gas condensing ou t of the hot p hase 
of the wind (iChevalier fe Cleggl 119851 : IWand Il995l) . In 
support of this scenario, we demonstrated that the lu¬ 
minosity of this component increases steeply with the 
SFR, following predictions for the total power radiated 
by the hot wind. While most of that energy would be 
radiated in X-rays, we belive we have identified the few 
times 0.1% of that cooling radiation emitted in Lya and 
optical lines (by its large Doppler shift). Since the cool¬ 
ing timescale will be decreasing rapidly as the wind tem¬ 
peratures falls to where the gas emits strongly in these 
lines, we expect it to be unstable and forms dense chimps 
or filaments. On qualitative grounds, some association 
with fast shocks seems likely; however it remains to be 
shown (using numerical simulations for example) why the 
shock speeds are lower than the fastest outflowing gas. 
Indeed, the low-ionization-state absorption is widely be¬ 
lieved to mark outflowing clouds accelerated by the hot 
wind and/or radiation; and the shocks could be gener¬ 
ated by the interaction of this low-ionization-state com¬ 
ponent with the cooling hot wind. 

If the broad emission-line wings do in fact mark the 
cooling hot wind, then the absence of this feature in the 
spectra of many starburst galaxies requires explanation. 
While fully exploring this topic is beyond the scope of 
this paper (see Thompson etal. , in prep), the winds 
from the galaxies with the highest SFR surface densities 
would appear to be good candidates for small cooling 
radii. Denser winds cool faster and reach the critical 
inflection in the cooling curve closer to the starburust 
region where the emission is brighter due to the higher 
density (prior to geomtrical dilution by the expansion). 

The relationship between the various outflow diagnos¬ 
tics in ULIRGs should be applicable to a wide range of 
intensely star-forming galaxies over cosmic time. These 
ULIRGs share a common property with high-redshift 
galaxies, namely high gas accretion rates, and evolve to¬ 
wards a quiescent state on the merger timescale. When 
ordered in a temporal sequence based on the merger 
progression, we find the fastest outflows and the low¬ 
est columns of neutral gas in the most advanced merger 
stages. These late stages show a larger contribution to 
Lfir from AGN, but the outflows develop long before 
the AGN dominates Lfir • The combination appears to 
drill holes through the neutral gas which allow a small 
percentage of the Lyman continuum to escape from even 
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these extremely dusty galaxies. 

One of the important outcomes of this study is that 
the prominent blue wings on the Lycc profiles also ex¬ 
ist on the optical emission lines; the optical wings are 
less prominent but cover the same velocity range as the 
Lya wings. It should be possible to resolve and study 
these fast outflows over much of cosmic time using highly- 
mutiplexed, infrared spectroscopy from the James Webb 
Space Telescope. Deeper observations with COS could 
provide detections of resonance absorption lines across a 
broad range of ionization states further defining the rela¬ 
tionship between the signatures of outflows in absorption 
and emission lines. 
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TABLE 1 

ULIRG Properties 


Name 

(i) 

z 

(2) 

log L fir 
(Lq) 

(3) 

F v (25[im)/ 

Fv(60f-im) 

(4) 

jpMW 
^B-V 

(5) 

rriNUV 
(AB mag) 
(6) 

m F uy 
(AB mag) 
(7) 

Merger 

Class 

(8) 

Separation 

(kpc) 

(9) 

IRAS 00262+4251 

0.097240 a 

11.90 

0.11 

0.066 

18.73 

19.22 

IVat 

< 1.4 

IRAS 01003-2238 

0.117701 b 

11.94 

0.25 

0.015 

18.49 

18.51 

vt 

< 3.4 

IRAS 08030+5243 

0.083500 3 

11.82 

0.07 

0.037 

21.78 

21.17 

IVbt 

< 1.3 

IRAS 09583+4714 

0.085887 c 

11.78 

0.18 

0.012 

18.46 

19.09 

Illat 

14.5 

IRAS 10378+1109 

0.136238 b 

12.08 

< 0.15 

0.028 

21.04 

21.22 

IVb 

< 1.2 

IRAS 11598-0112 

0.150694 b 

12.25 

< 0.22 

0.023 

17.54 

17.87 

IVb 

< 1.6 

IRAS 12071-0444 

0.128360 b 

12.09 

0.21 

0.039 

19.46 

20.82 

IVb 

2.25 

IRAS 16474+3430 

0.111471 b 

11.96 

0.11 

0.002 

18.39 

18.93 

Illb 

3.73 

IRAS 16487+5447 

0.103600 a 

11.98 

0.063 

0.013 

18.16 

18.47 

mb 

5.16 

IRAS 17028+5817 

0.106090 b 

11.99 

0.032 

0.020 

20.26 

20.82 

Ilia 

25.3 

IRAS 23365+3604 

0.064480 a 

11.96 

0.11 

0.097 

18.91 

19.88 

mb 

< 1.1 


Note. — Col 1: Name of COS target. Col 2: Redshift. Col 3: Far-infrared luminos¬ 
ity defined as L FIR = 3.86 x \Q b d 2 Mpc [2.58F u (60nm) + F v (100 iim)]L & . The fluxes (in Jan- 
skys) come from iMurphv et al.l (1199611 . and the luminosity distance assumes a universe with 
Oa = 0.7, f l m = 0.3, h = 0.7. In the limit that stars, rather than an AGN, heat the dust which emits 
the fa r-IR radiation, the inferred star formation rate is SFR = Lpj R /5.8 x 10 9 M 0 vr~ 1 (|Kennicuttl 
i 19981 ). The bolometric luminosities, Lj Ri are a factor of two larger than Lfir■ Col j: Infrared 
color. Ratio computed from th e flux densities at 25 /imand 60 /unlisted in Table 1 of lMurphv et al.l 
(|1996H . I Sanders et al.l (|1988bH used this ratio to identify ULIRGs with unusually high dust tem¬ 
perat ure; The warm ULIRGs hav e 0.2 < F v (25/tm)/ F v (60/tm) < 2.0. Col 5: Fore ground extinc- 
tion (Sch laflv fe Fi nkbcincr 201lt) reca librated from the infrared-based dust map of Schl eeel et all 
(|1998 1. The recalibration assumes a iFitzpatrickl (|1999f ) reddening law with Ry = 3.1. Col 6: 
Near-UV magnitude from GALEX (not corrected for foreground extinction). Col 1: Far-UV mag¬ 
nitude from GALEX (not co rrected fo r for egroun d extinction). Col 8: Interaction classification as 
described in Section 3.4 of lYeilleux et al.1 (120021) : I - First approach, II - First contact, III(a/b) 
- Pre-merger (Wide binary/Close bin ary), IV(a/ b) - Merger (Diffuse/Compact), V - Old Merger. 
Dagger ind icates object c l assifie d bv ISoto fe Martini (I2012H . Col 9: Separation of nuclei. Upper 
limits from IMurphv et abl (|1996f l: all others made from the observations presented here. 
a Redshift from CO in lSolomon et al.l (119971 ) 
b Redshift from optical emission lines. 
c Redshift from infrared Paa emission line. 
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TABLE 2 

COS Imaging Observations of ULIRGs 


Program 

Visit 

HST ID 

Config 

Filter 

Target 

Date 

T 

GO HST 







(s) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

13407 

1 

JC8T01010 

ACS/SBC 

F125LP 

IRAS 00262+4251 

2013-12-28 

892 

13407 

1 

IC8TA1010 

WFC3/UVIS 

F225W 

IRAS 00262+4251 

2013-12-28 

892 

13407 

1 

IC8TA1YMQ 

WFC3/UVIS 

F625W 

IRAS 00262+4251 

2013-12-28 

71 

5982 

n/a 

U2020106T 

WFPC2 

F439W 

IRAS 01003-2238 

1995-10-21 

800 

12533 

21 

JBRM21010 

ACS/SBC 

F125LP 

IRAS 09583+4714 

2012-06-08 

934 

12533 

21 

IBRMA1010 

WFC3/UVIS 

F225W 

IRAS 09583+4714 

2012-06-08 

704 

12533 

21 

IBRM1AM0Q 

WFC3/UVIS 

F625W 

IRAS 09583+4714 

2012-06-08 

75 

12533 

9 

JBRM09010 

ACS/SBC 

F125LP 

IRAS 12071-0444 

2012-05-17 

934 

12533 

9 

IBRMA9010 

WFC3/UVIS 

F225W 

IRAS 12071-0444 

2012-05-17 

704 

12533 

9 

IBRMA9F2Q 

WFC3/UVIS 

F625W 

IRAS 12071-0444 

2012-05-17 

75 

12533 

22 

JBRM22010 

ACS/SBC 

F125LP 

IRAS 16474+3430 

2012-12-14 

934 

12533 

22 

IBRMA2010 

WFC3/UVIS 

F225W 

IRAS 16474+3430 

2012-12-14 

710 

12533 

22 

IBRMA2F2Q 

WFC3/UVIS 

F625W 

IRAS 16474+3430 

2012-12-14 

75 

13407 

2 

JC8T02010 

ACS/SBC 

F125LP 

IRAS 16487+5447 

2014-01-12 

932 

13407 

2 

IC8TA2010 

WFC3/UVIS 

F225W 

IRAS 16487+5447 

2014-01-12 

932 

13407 

2 

IC8TA2NDQ 

WFC3/UVIS 

F625W 

IRAS 16487+5447 

2014-01-12 

91 

12533 

23 

JBRM23010 

ACS/SBC 

F125LP 

IRAS 23365+3604 

2012-09-23 

934 

12533 

23 

IBRMA3010 

WFC3/UVIS 

F225W 

IRAS 23365+3604 

2012-09-23 

704 

12533 

23 

IBRMA3PEQ 

WFC3/UVIS 

F625W 

IRAS 23365+3604 

2012-09-23 

76 


Note. — Col 1: HST Guest Observer program number. Col 2: Visit number. Col 3: 
Program identifier. Col 4 : Instrument configuration. Col 5: Imaging filter. Col 6: Target 
name. Col 7: Date Col 8: Exposure time. 



TABLE 3 

COS G130M Observations of ULIRGs 


Program 

Visit 

HST ID 

Target 

Date 

Acquisition 


RA 

DEC 

CWL 

Band (A/B) 

Duration 

GO HST 










(A) 

(A) 

(s) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


(7) 

(8) 

(9) 

(10) 

(11) 

12533 

1 

LBRM01010 

IRAS 00262+4251 

2011-11-11 

peak 

00 

28 

54.1526 

+43 08 15.82 

1291 

1137-1274/1292-1432 

592 

12533 

1 

LBRM01020 

IR AS 00262+4251 

2011-11-11 

peak 

00 

28 

54.1526 

+43 08 15.82 

1309 

1154-1293/1310-1448 

592 

12533 

2 

LBRM02010 

IRAS 01003-2238 

2011-12-03 

peak 

01 

02 

49.9631 

-22 21 57.02 

1291 

1137-1274/1292-1432 

1716 

12533 

4 

LBRM04010 

IRAS 08030+5243 

2011-10-30 

peak 

08 

06 

50.8856 

+52 35 7.70 

1291 

1137-1274/1292-1432 

932 

12533 

18 

LBRM18010 

IRAS 09583+4714-B 

2012-11-08 

peak 

10 

01 

30.4490 

+46 59 51.74 

1291 

1137-1274/1292-1432 

1832 

12533 

6 

LBRM06010 

IRAS 10378+1109 

2011-11-11 

peak 

10 

40 

29.1743 

+10 53 18.17 

1291 

1137-1274/1292-1432 

420 

12533 

6 

LBRM06020 

IRAS 10378+1109 

2011-11-11 

peak 

10 

40 

29.1743 

+10 53 18.17 

1309 

1154-1293/1310-1448 

420 

12533 

8 

LBRM08010 

IRAS 11598-0112 

2011-11-18 

peak 

12 

02 

26.7505 

-01 29 15.49 

1309 

1154-1293/1310-1448 

652 

12533 

8 

LBRM08020 

IRAS 11598-0112 

2011-11-18 

peak 

12 

02 

26.7505 

-01 29 15.49 

1327 

1173-1312/1328-1468 

652 

12533 

17 

LBRM17010 

IRAS 12071-0444 

2013-01-24 

offset 

12 

09 

45.1000 

-05 01 13.20 

1291 

1137-1274/1292-1432 

588 

12533 

17 

LBRM 17020 

IRAS 12071-0444 

2013-01-24 

offset 

12 

09 

45.1000 

-05 01 13.20 

1309 

1154-1293/1310-1448 

588 

12533 

24 

LBRM24010 

IRAS 16474+3430 

2013-04-14 

offset 

16 

49 

14.2380 

+34 25 8.68 

1291 

1137-1274/1292-1432 

2204 

12533 

13 

LBRM13010 

IRAS 16487+5447 

2011-11-06 

peak 

16 

49 

47.0306 

+54 42 35.41 

1291 

1137-1274/1292-1432 

1888 

12533 

14 

LBRM14010 

IRAS 17028+5817 

2011-10-26 

peak 

17 

03 

42.0084 

+58 13 44.39 

1291 

1137-1274/1292-1432 

1456 

12533 

16 

LBRM16010 

IR AS 23365+3604 

2011-11-07 

peak 

23 

39 

1.2726 

+36 21 8.45 

1327 

1173-1312/1328-1468 

1340 

12533 

20 

LBRM20010 

IR AS 23365+3604 

2012-11-24 

offset 

23 

39 

1.2290 

+36 21 8.03 

1327 

1173-1312/1328-1468 

1876 


Note. — Col 1: HST Program ID. Col 2: Visit number. Col 3: HST ID. Col Target name. Col 5: Date of observation Col 6: Type of acquisition: 
blind offset or peak up. Col 7,8: Coordinates. For blind offset acquisitions, these coordinates are the exact position of the PSA. For peak ups, the telescope was 
pointed near these coordinates and then moved in a spiral pattern to cover a square grid (ACQ/SEARCH with SCAN-SIZE = 2 and STEP-SIZE = l") recording 
a spectrum at each stop. The telescope was then moved to the flux-weighted mean position of the target. We then improved the centering in the cross-dispersion 
direction by taking a spectrum, collapsing it in the dispersion direction, finding the offset in the cross-dispersion direction to the mean of the collapsed spectrum, 
and finally slewing the telescope by this offset. After this ACQ/PEAKXD, we moved the telescope along the dispersion direction to obtain several ACQ/PEAKD 
exposures (STEP-SIZE=1.118; CENTER FLUX WEIGHTING) and found the center of the knot along the dispersion direction. The science spectra were then 
obtained through the PSA. Col 9: Central wavelength. Col 10: Wavelength coverage. Each spectrum has a gap between detector Segments B and A. For some 
targets, we simply chose a CWL that avoid loosing important lines to this gap, but for others we needed continuous coverage and requested exposures at two 
CWL’s as indicated. Col 11: Exposure time. 
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TABLE 4 

Measured Properties from COS Spectra 



Object 

F(Lya) x 10 15 

logL(Lya) 

/esc(Lya) 

W r (Lya) 

E(B-V)* 

Vs 20 

V%20 

W r (Si II) 

14+(Si III) 



( ergs s —1 cm -2 ) 

(ergs s -1 ) 


(A) 

(mag) 

( km s _1 ) 

( km s -1 ) 

(A) 

(A) 


(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

IRAS 

00262+4251 

q cq0.53 

O.OZ — o 61 

40.92 

2.05lo34 x 10- 2 

—10.8±|'| 

O qo+0.03 
u.oz_ 0 03 

-326±4 3 

449l 3 g 

1.04 

0.93 

IRAS 

01003-2238 

70 116.6 
' y - i -4.00 

42.53 

7.2 1° 0 \ x 10" 2 

-32.4l 45 0 

0 02"*" 0 02 
U - [JZ ~0.02 

-1002l 7 ® 8 

692 I 4 Q 4 



IRAS 

08030+5243 

< 0.072 

< 39.09 

< 1.2 x 10' 4 







IRAS 

09583+4714-B 

2-91°jf 8 32 

40.73 

i.3j;°- 2 x 10- 2 

7 q+3.6 
' '^—4.0 


91 9+48 
Z 1 Z _ 27 

500± 2 



IRAS 

10378+1109 

c oo 0.12 
o.oo—o 92 

41.42 

2 - 44+ n2i 6B x 10-2 

< -13 


-664± 275 

440j; 2 2 



IRAS 

11598-0112 

151 121 

±0 -33.0 

< 43.00 

5 5 +4 - 4 
o.o_i 2 

-30.itjg;5 

oi 

—6679” 1 " 2137 
oo/y_2408 

4291^203 


0.05 

IRAS 

12071-0444 

5 qq 0.88 

1.03 

41.41 

0.241^1 

< -17 


_743+ 153 

ffto -120 

33 it 28 



IRAS 

16474+3430 

2 31 101 

40.87 

l-33±g;gg X 10- 2 

— 1 2 +1 ' 7 
1,z -3.4 

0 16+ 0 - 02 
u,iD —0.13 

-97t?oi 

97 +I 8 

-27 

1.32 

3.42 

IRAS 

16487+5447 

0 7«1.86 
z -'°_0. 87 

40.88 

i-3i±8;i| X 10 - 2 

9 o+1.5 
z.Z —2 1 

0 12 - *" 003 
1 ^_0.02 

—408^427 

744±t 2 

0.91 

1.01 

IRAS 

17028+5817 

< 0.059 

< 39.23 

< 4.6 x 10 ~ 4 







IRAS 

23365+3604 

5-69™ 7 75 

40.78 

2.10l“ ° 8 x 10~ 2 

< -29 


-1046l^ 7 

146t^ 




Note. — Col 1: Name Col 2: Lya flux. The uncertainties include a term for the continuum placement. We adopted the linear continuum 
fit to bandpasses on both sides of the line and compared this to the maximum and minimum values of the continuum one could fit by adopting 
different bandpasses, e.g., such as our global fit to the continuum. Col 3: Lyaescape fraction. Col 4 : Lyaluminosity Col 5: Equivalent width 
of Lya. Where the continuum was not detected, we use the 3<r upper limit on the continuum, as measured in a 20A wide window, to estimate 
the minimum equivalent width (which appears as an upper limit since negative W r indicates net emission). Col 6: Stellar reddening. Values 
assume an intrinsic far-UV continuum slope of /3 = —2.3 ± 0.1, where F\ oc A^ and the starburst reddening curve of ICalzetti et al.l (1200011 . 
The continuum bandpasses regions were chosen to avoid telluric emission, foreground Galactic absorption lines, and interstellar absorption from 
resonance transitions. Col 7: Interpercentile velocity at p = 20% for the blue wing (relative to the systemic velocity. Col 8: Interpercentile 
velocity at p = 20% for the red wing (relative to the systemic velocity. The line width is IP \20 = Vb 20 + Vr 20 - Col 9: Absorption equivalent 
width of Si II 1190. The stronger Si II 1260 line is often blended with S II 1259. Col 10: Absorption equivalent width of Si III 1206. 













TABLE 5 

Measured Properties from Optical Spectra 


Object 

E(B-V) 

F(Ha) x 10 15 

Aper. 

11+(Na I) 

Vb2o(5007) 

Vr2o(5007) 

Vsh 


(mag) 

( ergs s —1 cm -2 ) 

Corr. 

(A) 

( km s -1 ) 

( km s -1 ) 

( km s -1 ) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

IRAS 09583+4714-B 

0.138 + 0.003 

4.63 + 0.01 

1.81 

0 

-48t“ 

72l° 

£ 150 

IRAS 16487+5447 

0.115 + 0.001 

13.75 + 0.01 

1.31 

< 0.5+ 

—217t?j 

238 + o 2 

220 ± 10 

IRAS 16474+3430 

0.199 + 0.001 

7.849 ± 0.005 

1.51 

2.53 

-80t? o 

35llJ° 

300 ± 25 

IRAS 00262+4251 

0.242 + 0.004 

6.466 + 0.007 

1.51 

1.0 ±0.3+ 

—5031?? 

1721^ 

225 ± 50 

IRAS 23365+3604 

0.318 + 0.006 

7.91 + 0.02 

1.55 

4.5 ±0.4+ 

—687 tir 

1501? 

150 ± 50 

IRAS 10378+1109 

0.339 ± 0.002 

6.600 ± 0.004 

1.64 

6.7 

-671±|| g 

2541? 

320 ± 10 

IRAS 12071-0444 

0.098 + 0.001 

1.578 + 0.001 

1.33 

1.45 

-7231?! 

1641J 1 

350 ± 20 

IRAS 01003-2238 

0.047 + 0.004 

6.73 ± 0.02 

1.64 

0.77 

-16501J 0 

68l? 

260 ± 10 

IRAS 11598-0112 


1.76 + 0.02 

1.71 

0.53 

—18261® 4 ® 

6521?? 

300 ± 15 

IRAS 08030+5243 

0.513 + 0.002 

8.700 ± 0.005 

2.27 

7.4+ 1.0+ 



150 ± 100 

IRAS 17028+5817 

0.143 + 0.007 

6.485 ± 0.008 

1.56 




225 ± 15 


Note. — Col 1: Target name. Col 2: Nebular reddening inferred from Balmer decrement measured in 
ESI spectrum. We adopted an intrinsic flux ratio F(Ra)/F(H/3) = 2.86 appropriate to Case B recombination 
at an electron temperarture of T e = 10 4 K (IQsterbrock fc Ferlandll2006h and used the ICardelli et aTI (119891) 
reddening curve with a ratio of total-to-selective extinction Ry = 3.1. We extracted the ESI spectrum from 
2 , . , 5 long aperture centered on the COS G130M pointing. The optical emissio n lines were fit si multaneously 
using two emission components plus stellar Balmer absorption as described in ISoto et al.l (120121 ). Col 3: Ha 
flux in 2 , /5 wide apeture along the ESI slit. Col 4'- Aperture correction required to match Haaperture to COS 
G130M Lyaaper ture. We used a n r-band image of each ULIRG as a template for the spatial distribution of 
the Haemission (IKennicuttll 19891) . We scaled the relative intensities by the COS vignetting function described 
in the STSci Instrument Science Report COS 2010-09(vl); this correction lowers the weight of flux at radii 
greater than approximately O'/T We define the aperture correction by the ratio of the vignetted flux in the 
COS aperture (2'/5 diameter circle) to the ESI apert ure (2 , /5 by 1" 0 rectangle). Col 5: Absorption equivalent 
width of Na I. A dagger denotes measurements from IMartinl (120051 ); all other measurements were made from 
the ESI spectra using the method described therein. Col 6: Interpercentile p = 20% velocity (blueshifted 
relative to systemic). Col 7; Interpercentile p = 20% redshifted velocity (redshif ted relative t o syst emic). Col 
8: Shock velocity fitted to F([OIII])/YL/3 vs. F([OI ])/Ha flux ratio using the lAllen et alj (|2008l ) library of 
fully radia tive shoc k mode l s calc ulated with the MAPPINGS III shock and photoionization code. We list the 
value from ISoto fc Martini (120121 ) aperture that falls within the COS primary science aperture. 
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TABLE 6 

Radiative Transfer Models for a Neutral Shell 


Name 

log N hi 
( cm -2 ) 

V 

(km s -1 ) 

T~d 

fc 

IRAS 00262+4251 

19.9 

0 

0 

0.3 

IRAS 01003-2238 

17.0 

200 

0.01 

0.5 

IRAS 09583+4714B 

20.0 

-80 

0.4 

0.9 

IRAS 10378+1109 

19.5 

-30 

0.01 

0.5 

IRAS 11598-0112 

18.8 

±150 

2.8 

1.0 

IRAS 12071-0444 

20.3 

0 

0 

0.8 

IRAS 16474+3430 

18.5 

+250 

2.0 

1.0 

IRAS 16487+5447 

19.5 

-50 

0.6 

1.0 

IRAS 23365+3604 

19.6 

-400 

0.05 

1.0 


Note. — Col 1: Name. Col 2: Column den¬ 


sity of neutral hydrogen. Col 3: Shell velocity 
using the convention that outflows have a neg¬ 
ative sign. Col 4 : Dust optical depth. Col 5: 
Covering fraction of neutral gas. 
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APPENDIX 


VELOCITY MEASUREMENTS 

In order to compare Lya profile to the optical emission lines in this work, we require q uanti t ative measurements of 
the shapes of the emission-line profiles. Below, we follow the method introduced by i Whittle! (|1985f) to describe the 
[O III] profiles of active galaxies. This quantitative descriptions should also facilitate comparison of the highly unusual 
Lya profiles of ULIRGs to those of other samples of Lya emitting galaxies. 

The strategy rests on measuring line widths at a specified fraction of the total profile area. We take the profile 
area to be the net area in emission where the integration is over the entire line profile; any other definition would be 
sensitive to spectral resolution. For purposes of measuring the line width, asymmetry, etc., one chooses a value of the 
fractional area that isolates the desired profile feature such as the line core, wings, etc. 

We call the width at a specific fractional profile area the interpercentile velocity (IPV). In summary, IPX 10 , IPX' 20 , 
IPX' 30 , •••., and IPV 90 mark the bandpasses that include, respectively, 90%, 80%, 70%, ..., and 10% of the profile 
area. To measure IPV. 20 , for example, we integrate the velocity profile up to a velocity V B that includes 10% of the 
total area and then continue integrating up to a velocity, V Rl that includes 90% of the total area. The interpercentile 
velocity width IPV p includes the central 1 — p percent of the line flux, where p de note s t he f raction the line emission 
detected at velocities larger than the interpercentile velocity width. Following IWhittlei (Il985f) . we define 


and an asymmetry parameter 


IPV p = v Bp + v Rp 


_ v Bp — v Rp 
Lp - v Bp + V Rp 


(Al) 

(A2) 


When applying this scheme to Lya, or any line that have absorption as well as emission components, we have chosen 
to define the area of the line profile as the total area such that absorption components subtract from the total area. 
We estimate the uncertainty in each marker by varying the continuum placement as described in Section [2] 

In Figures a-1, we have marked velocities of the blue wing, V Bp , and red wing, V Rp , where 20%, 50%, and 80% of the 
total emission in the line lies at higher velocity. These two markers are defined for any percentile p of the total profile 
area, so we show how the line width grows as p decreases to include more of the total profile area. The Doppler shifts 
of the blue and red wings are generally not symmetric about either the systemic redshift or the median velocity the 
profile area. When the profile is entirely in absorption blueward of the systemic velocity, as for example in Figure c, 
then the velocity Vr 20 will be redward of the systemic velocity because it is defined in terms of the net profile area in 
emission. 

Since we have defined profile asymmetry relative to the median velocity of each emission-line profile, a profile with 
a net Doppler shift can have A ss 0 (no asymmetry) if the two sides of that profile are symmetric. Returning to the 
P Cygni profile in Figure c as an example, again, we see the emisssion is only mildly asymmetric about the median 
velocity of 300 km s , and the Lya asymmetry parameter is driven slightly redward by the presence of the blueshifted 
absorption. This gross asymmetry with respect to the systemic redshift is recognized instead by the comparison of 
V BtP and Vr iP which have the same sign instead of opposite signs. 

As a guide towards interpreting these measurements, we provide a few brief comments on the individual objects. 


IRAS00262+4251 — Black Lya absorption near the center of the [O III] emission profile. Scattered Resonance photons 
broaden Lya. Despite the disimilar line shapes, the blue marker velocities have similar values in both lines. The [O III] 
profile shows a distinct blue asymmetry in the line wings whereas the Lya profile has a red asymmetry at percentile 
markers. 


IRAS 01003-2238 — This object has a substantial AGN luminosity. The Lya profile has a blueshited (-200 km s *) 
median velocity about which it is fairly symmetric. The [O III] emission has a net blue asymmetry in the line wings. 

1RAS09583+4714-B - This non-ULIRG spectrum shows the classical P Cygni profile as described in the text. The 
asymmetry parameter of the Lya profile is not very strong because it is measured relative to the median velocity which 
is redshifted 300 km s -1 . 


IRAS 10378+1109 — The blue wings of the Lya and [O III] profiles share the identical shape blueward of -200 km s 1 . 
The Lya profile is missing photons in its core relative to the [O III] profile but shows more redshifted emission. 

IRAS 11598-0112 — We show this narrow-line Seyfert 1 whose UV spectrum has mini-BAL characteristics for com¬ 
pleteness. The [O III] profile blends with Fe II lines, so we compare to H/3 instead. The H/? profile shows the blue 
asymmetry that we typically find on optical lines from ULIRGs. The large Lya velocities are probably misleading 
because Lya is likely blended with metal lines. In the absence of blending, this spectrum would be the only one with 
Lya emission extending to larger blueshifts than the optical lines. 
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IRAS12071-04H — Resonance photons scatterd out of the deep absorption trough 200 km s -1 emerge at bluer and 
redder wavelengths. The bluesliifted component has a prominet wing which also shows optical emission. The Lya and 
optical line shapes are very similar blueward of -500 km s -1 . 

\RAS16i7^+3430 — The COS aperture includes a second nucleus offset 1"84 from the target. In the ESI spectrum, this 
offset nucleus is redshifted about 300 km s _1 relative to the target. Although the COS spectrum contains no spatial 
information, the redshifted Lya absorption trough may be associated with the offset nucleus. The fast orbital motion 
of these nuclei explains why the optical emission lines are broad and double peaked. The lack of Lya emission from 
the offset nucleus explains why the Lya profile is narrower than the optical lines. 

IRAS 16487+5447— Much of the Lya profile shows absorption. The emission is largely redshifted. With a median 
velocity of roughly 500 km s -1 for the net emission, the weak, the blueshifted emission at -300 to -500 km s -1 
produces a strong blue asymmetry as shown by the IPV 50 and IPV 20 markers. 

IRAS 23365+3604 — The Lya profile is not delineated by narrow absorption components, yet the larger amplitude of 
the peak [O III] emission suggests some scattering out of the line core. The blue wings of the Lya and [O III] lines 
cover the same velocity. The extra wiggles in the blue wing of the Lya result from absorption by an expanding shell 
in the radiative transfer model. 



35 



0.5 

< 0 
-0.5 


‘"lOOO 

a 

" 500 

> 

P-, 

” 0 

0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 

Fraction of Profile Area at |V| > IPV Fraction of Profile Area at |V| > IPV 

(a) IRAS00262+4251 



500 


> 

>-500 


& 

A 


^ 1 

^ 4 

A 4 


8 1 


* * A 

A * * A 

s ■ 1 * 


Fig. 1.— Top Row: The Lya and [O III] line profiles normalized by the net flux in emission. For the Lya line, we have marked the 
interpercentile velocities at p = 20% (magenta), p = 50% (yellow), and p = 80% (cyan) to illustrate the meaning of Vb, p and Vr, p - Bottom 
Row Left: Interpercentile velocity width and asymmetry parameter vs. fractional area of the profile for Lya and [O III]. Note that the 
asymmetry is measured relative to the median velocity which is marked by the dashed line in the top and right panels. Bottom Row, Right: 
Comparison of the velocity markers Vb, p and Vr, p for Lya and [O III]. The blue wings (squares) are more correlated than the red wings 
(triangles). Blue and white points denote [O III] and Lya, respectively; and the open symbols denote upper and lower limits on the solid 
point. 
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(b) IRAS01003-2238 
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(c) IRAS09583+4714-B 
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(d) IRAS10378+1109 
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(e) IRAS11598-0112 
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(f) IRAS 12071-0444 
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(g) IRAS16474+3430 
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